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ABSTRACT 


Lumped-parameter and finite-element-models are devel- 
oped for the dynamic behavior of plane straight sided 
frames. The models do not include axial and shear defor- 
mation, but the equations of motion developed using the 
models allow for time varying external loading. The per- 
formance of these models is evaluated by a comparison with 
a standard transfer matrix method for the special case of 
free undamped vibration. . The Figite-element-model proves 
to be much the better model. For the first five modes, 
the finite-element-model with 27 degrees of freedom differs 
by no more than 0.4% from the "exact values" given by the 
transfer matrix method. The lLumped-parameter-model with 
39 degrees of freedom gives errors roughly ten times as 
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NOMENCLAUTRE 


Cross sectional area of member i. 
Spectral matrix from diagonalization of J. 
Matrix relatwnge@d co Cc. 

Coefficient vector. 

Element of coefficient vector 

Modal matrix from diagonalization of J. 
Deflection and slope vector. 

Modulus of elasticity of member i. 
Flexural rigidity 

External force vector. 

Vector of x components of external force. 
x component of external force at node j. 
Vector of y components of external force. 


y component of external force at node j. 


Row vector in powers of r for —=. 


Row vector in powers of r for $, 


Identity matrix. 


Second moment of area A; with respect to a centroidal 


axis perpendicular to plane of frame. 
Matrix sum of inertia terms. 


Inertia matrix. 


Mass moment of inertia of segment j about an axis 
normal to the plane of the structure and thru the 


segment centroid. 
Rotational inertia between j and jtl. 


Stiffness matrix. 


K 


Is 


Stiffness of segment j. 


Incidence matrix. 


Ss 
Length of segment of member i, = — P 


Mass matrix. 

Rotational mass matrix. 

Translational mass matrix. 

segment mass. 

Moment vector. 

Moment applied between node j and j-l. 
Total number of subdivisions. 

Number of subdivisions in member i. 


Matrix used to apply end conditions in lumped- 
parameter-model. 


Row matrix in powers of r for w. 

Moment arm matrix for external force vector. 
Generalized force vector. 

Corner transformation matrix. 


Axial coordinate. 


Matrix@prcauce of 5 
Length of member 1. 
Kinetic energy. 
Rotational kinetic energy. 
Translational kinetic energy. 
Potential energy. 

Vector of x components of displacement 


Axial displacement. 
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us = x displacement of node j. 
Vv = Vector of y components of displacement, 
vs = y displacement of node j. 
Ww. = Specific weight of member i. 
Ww = Transverse deflection. 
x = Matrix of the x components of distance. 
. = x-coordinate of the centroid of segment j. 
ij = X component of distance from node j to node it+l. 
¥ = Matrix of the y components of distance. 
us = y-coordinate of the centroid of segment j. 
ae = y component of distance from node j to node itl. 
1 
Zz = product of ae and D* - 
(x,y,2) = Cartesian coordinate system. 
O. = Deviation angle at a corner of the frame. 
s) = cos Qa. 
Y = Angle between axial direction in local coordinate 
system and x direction in Cartesian system. 
n = sin a. 
0 = Relative rotation vector. 
F = Relative rotation of hinge j. 
Ul = Mass per unit length. 
> = Slope. 
Ww) = Natural circular frequency. 
Subscripts 
a = assembled matrices. 
a = Property of member 1. 
5 = Property at node j. 


deli 


ij = Property from node i to node j in finite-element- 
model; from node j to node i+l in lumped-parameter- 


model. 

u = Upper rows of a matrix as determined by end con- 
Gikeaoniss 

i = Remaining rows of matrix. 


Speers eElpiEs 


te = Transpose of vectors and matrices. 


: = Matrix after application of end conditions. 
: = Vector resulting from corner transformation. 
= In lumped-parameter-model denotes a once reduced 


matrix; in finite-element-model denotes the proper- 
ties of the element to the right of the corner. 


le 


INTRODUCTION 


When designing a structure, the engineer must be aware 
of the many dynamic disturbances that may act upon this 
structure. He must also be aware of the response that the 
structure will have to these disturbances. The calculation 
of the response is usually a very difficult undertaking. 

To assist in this calculation, the structure is generally 
characterized by a simplified model, -and the analysis 
Carried out uSing this model. 

In this paper, two models for analyzing plane frame 
structures with "chain" topology* are developed. The 
first model lumps both mass and flexibility and the second 
is a finite-element-model. To test these models, analyses 
to find the in-plane circular frequencies and translational 
displacements are carried out on a typical gable bent. The 
results of these analyses are compared to frequencies and 
displacements obtained through the use of the so-called 
transfer matrix neehoasee 

The equations of motion as developed for the models 
provide for time varying external loads, but do not include 
damping. The analysis of the gable bent is carried out for 


the special case of free vibration. This simplification 


*Plane framed structures having 2 ends and a single 
path from end to end. 


**Superscripts refer to the bibliography which appears 
on page 48, 
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is required so that the principal modes obtained through 
the use of the models may be compared with values from the 


transfer matrix method. 
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DEVELOPMENT OF THE LUMPED—-PARAMETER-MODEL 


In order to find in-plane frequencies and displace- 
ments of a structure consisting of straight sides with uni- 
form cross section, a model lumping both mass and flexibil- 
ity is developed. This model neglects axial and shear de- 
formation, but can include rotatory inertia. 

The model is formed by subdividing each straight mem- 
ber into segments. These segments are considered to have 
square ends. The length along the centerline of the 
straight member and the sum of the lengths of the segments 
replacing it are the same, even though this would cause some 
overlapping at the corners. The mass and flexibility of 
each segment is lumped at the segment centroid. The struc- 
ture is then replaced by a chain of straight rigid bars, 
pin connected at the locations of the segment centroids. 
The mass of the segment remains concentrated at the hinge, 
and the flexibility of the segment is replaced by a tor- 
sional spring. The equations of motion are developed and 
the model analyzed to obtain its natural frequencies and 
the displacements of the hinges. 

Rigid body translation is prevented by taking the left 
end as fixed or pinned. The right end is allowed to assume 
any of the common end conditions: fixed, pinned, simply 
supported, or free. 

Development of the Equations of Motion 
To illustrate the method, a structure containing two 


straight members intersecting at an acute angle is used. 


Je 





(a) Two Member Frame 





a oe 


————_ 
(c) Lumped Model of 
Frame 


Figure 1 


Both ends of the structure are taken as fixed (Figure la). 


For each straight member i (i=1,2), let 


= 
il 


specific weight, lb/cu.in. 


j 

A, = cross sectional area, sq.in. 

I. = second moment of the cross sectional area with 
respect to the z axis, in.” 

EB. = modulus of elasticities, loys quae 

Ss. = length of member, in. 

n, = number of segments; nj =2, nj=3. 
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Each section ®f the structure “We dividee ito n, seg- 
ments (Figure lb). The n. need not be the same for the dif- 
ferent members, but in each member the segments are of uniform 
length. The parameters of each segment are then lumped at its 


centroid. These parameters are 


A = segment length = s./n; inches. 
Tics segment mass = (AWL) » /386 ecchs/inn 
J; = mass moment of inertia abhout an axis normal 


to the plane of the structure and through the 
segment centroid, 


2 ; 2 
m. (2/12 + I,/A;) in-lb-sec. 


k, = segment stiffness = (EI) -/2; lb-in. 


The structure is then replaced by the rigid bars and tor- 
sional springs (Figure lc). The springs are unstrained in the 
datum position. The Cartesian coordinates of the end points 
and hinges of the structure in the datum position are taken 
as 2a and Y4 where j=1,2,3,4,5,6,7. Thus the displacements 
of the hinges ee: Ya in the positive x and y directions for 


a relative rotation, 8., at hinge j are given by 


J 
us -¥ @ 
cr) 
ia | 
where 
ue col(u,, Uz, Ug, Up, Up, us) 
Vv = col(v,, V3r Var Vor Ver V5) 
68 = col(®,, 85, O., Oy, On, 86) 


i 


a me years 
Yo] Y59 0 0 0 0 
J esme pre ecy gat eens 
ye 
= Yaa Yap, Yate. 
( 
ia Y52 Yoo Ygumisgal 
Yer Ye2 Ye3 Yea Yes 66 
a 
X54 X55 0 0 0 0 
X31 X45 X33 0 0 0 | 
x - | 
sid 4 
“41 “42 *43 “44 So 
} 
{ 
x50 *Somm5 3 ~5a8 hae? 
{ 
i *61 NG 63 764 “65 "66 
oe isl ey 
i Oa i 


The following notation has been used in this equation: 

1. Column vectors are denoted by lower case symbols. 

2. Rectangular matrices are denoted by capitalized symbols. 

3. All vectors and matrices are underlined. 

This notation will be used in all other such equations. 
Associated with each bar is a rotational inertia made up 

of half the mass moments of inertia of the preceding and fol- 


lowing segments. This rotational inertia is represented by 
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—_ 
c= x (J. + al: 
Using equation (1), the kinetic energy, T, and the poten- 


tial energy, U, of the structure are formed. 


. = =0°X0Mme + =—6 Yume S6CT a mo 
a cee ili eee a gg le 
Coy 
__ lg 
U = 59 KO 
where 
M = diag(m,, M3, My, Me, Me, 0) 
Jo= AiaglT,)+ Joo, Ja3r Jgge Jgg¢ Jeg? 
K = diag tig, ReRMe mc. , “hed 
i & eo ee me 0 | 
ali 1 0 0 0 7 
ul 1 A 0 0 0 
— 
1 1 ili i) 0 0 
if 1 a " 1 0 
i a} 1 a il. ii 


Derivatives with respect to time are denoted by dots and the 
transpose of a matrix by a superscript t. 


The Lagrangian equations of motion, 


May be written 


Jo + KO = q (3) 
where 
j= kK + Poy + Doo 
—— =o —_— —— —_—_- -—O-— 
eS = generalized force vector. 
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Let the x components of the external forces applied at hinges 
2 to 7 be represented by the vector 
ft, = col(f,.,, f 


f£ 1 


oo la ‘es ET mdse 
In a like manner define a vector io For the counterclockwise 
external moments ms applied between hinge i-1 and i, the vec- 


GOr 
n= col(m, m3) Mae Nese Me: m7) 
is defined. The generalized force vector may be expressed in 


terms of these vectors as 


Meee ot t 
g=Lnm ca Ee (4) 


it shouldsbe noted thatsthe matrices Quand K are bork 
symmetric. This condition results when any symmetric matrix 
is postmultiplied by a matrix and is premultiplied by its 
transpose. This symmetry property is highly desirable when 
dealing with matrix eigenvalue problems and should be preserved 
if possible when applying boundary conditions. 
Application of Boundary Conditions 

Up to this point the development has been general in na- 
ture. We now restrict further development to free vibration. 
With this restriction, the only external forces applied are 
those required to satisfy end constraints. If the model has 
both ends fixed, the end conditions and corresponding exter- 


nal forces are 


@D 
| 


i= 0, Mo x 0 
V4 = OF £05 # 0 
= 0, f7 # 0 


Oana o 


G 
~ 
| 


ee Ol cee 


Zz 5 4 


Equation (4) can therefore be reduced to 


qeQet (5) 
where 
£ = coliag £50 "hy, mt 
lent nye 
[0 X62 Yeo 1 
lO x y p 
| 3 ed 
Q = ; 
0 X64 Yea 4% 
O c= es me 


—_ 


jo “66 66 


The columns of Q are the first and” last columns of cae and 


= and der 


the last columns of X 
The end conditions may be represented by 


ove = 0. (6) 


This equation is derived from the following set of equations. 


G6, = 0 

uy = X61 974% 6999 t% 6 39 3t% 49 gtk en 9etX 0 = 0 
We. = Vig "1 goo ae ue ESS Vee Gl 
0, + 0, + 0, + 0, + 0, + O, = 0 


A process must be devised to incorporate the boundary 
conditions. This process will reduce the number of degrees of 
freedom by an amount equal to the number of constraint condi- 
tions, in this case four. The result will be a pair of homo- 


geneous equations with two remaining o are 


aa 


A normal Gaussian elimination technique for the solution 
of simultaneous equations is used to solve for fas ae oe 
and m- Through the use of the end conditions, the number of 
unknown Je is next reduced to two. The four non-homogeneous 
equations containing the forces and moments can be deleted, 
leaving only the two homogeneous equations. A similar elim- 
ination could be used for any other set of end conditions. 


This process can be represented by a set of matrix opera- 


tions. Equation (3) is written in partitioned form as 


| Fh 


(7) 


where the subscript u represents the rows which will contain 
the non-homogeneous equations, and & the remaining rows. This 
may require a row interchange to insure a non-singular O° 
In the present case om is 4x4. Gaussian elimination is now 
performed to solve for the forces and moments. In this proc- 


ess the elements of Q) will be reduced to zero. The lower 


portion of equation (7) may be rewritten 


J'6, + Ki6, = 0 (8) 
where the prime denotes reduced matrices. 


The end conditions, equation (6), may be rewritten 
a 
lot of] fo, 


and from this we get 


and therefore 
Lae, Beeebeat 


Thus ee may be eliminated from 6 by the transformation 


—2 = (9) 


{t+ 


where I is the identity matrix, in this case of order 2x2. 


It is easily shown that 


PQ = 0 
and that in equation (8) 
ic 
Nie coe 
Jo = Pg 
and 
K' = PK. 


Tes Oemay*thenebe répitacedmammequatiton™ (8) “by PO): After dif- 
ferentiating equation (9) twice, the 6 of equation (8) may 


also be replaced, and the following expression obtained: 


8, + Ki8, = 0 (10) 


en eas 
where 
a & 
Jo a P JP 
a t 
Ky 7 P KP 


The boundary conditions can therefore be applied and a set of 
homogeneous equations obtained by this process. This trans- 


formation does not destroy the symmetry of the matrices. 


Zz 


Eigenvalue Solution of the Equations of Motion 


Assuming that the vibration will be sinusoidal, equation 
(10) may be written 


it == Z te 
KpO, = w dhe, (an) 


where w is the natural circular frequency. This equation is 
in a form suitable for eigenvalue solution. 

The solution of equation (11) is now undertaken. This 
solution will involve eigenvalues and eigenvectors.” A sym- 
metric matrix may be expressed in terms of its modal matrix 
and its spectral matrix where the spectral matrix is a diag- 
onal matrix of eigenvalues and the columns of the modal matrix 
are eigenvectors. This property is used several times in the 
solution. A further important property of the modal matrix 
is that the eigenvectors are normalized with respect to the 
identity matrix. Thus the transpose of the modal matrix is 
also its inverse. 

af Jo is resolved into its (diagonal) spectral matrix 


Boand modal umaiesi1x D according, tomehe relation 


Up (21) ay 
equation (11) becomes 
- eZ 1c; 
Kio, = w DBD, . 
Since DD* = I, the above equation may be rewritten 


" t ao 2 |e 
me ey. 
Since B is a diagonal matrix of positive elements, we define 


1 
the principal square root as B2 where the elements of this ma- 


trix are the positive square roots of the elements of B. Not- 


ee 
ing that B 2B2 = I, the following form results, 
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—X = Q — Q 
= .t 
Le fbings Z=eB2p ci) , 
25 I 
Ky DB 22 = w DB2z 
2 
Solving for wz, 
i ae 
B2D°KsDB 22 = wiz . 
-i ft =8 
Letting S = B 2D KDB 2, the equation may be rewritten 
Sz = wz . 


By diagonalization of S, the values of ig and z may be ob- 
tained. From these values, the natural frequencies are ob- 
tained directly. The values of the rotations of the nodes 


1 
are obtained from the equation 8, = DB2z. The rotations 


Q 
at all the hinges of the structure can be recovered by use 
of equation (9). The translation displacements of the hinges 


are easily obtained with the use of equations (1). 
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DEVELOPMENT OF FINITE-ELEMENT-MODEL 
General 

A finite-element-model is developed to describe the 
in-plane motion of a frame with straight members of uniform 
cross section. The model neglects axial and shear deforma- 
tion and can include rotatory inertia if desired. 

The use of a finite element approach requires the divi- 
sion of the structure into elements, and the determination 
of element mass and stiffness matrices. The determination 
of the elemental mass and stiffness matrices will be taken 
up first. 

Consider a general element of length &, in a local 
coordinate system. In this system, r is the axial coor- 
dinate, w is the transverse deflection, and ¢ is the slope. 
The element has no flexibility in the r direction. If we 
assume that w may be approximated by a third degree polyno- 


mial in xr, then 
3 


HN 


2 
WwW Co + cir + Cor + Cor 
and 


dw _ 2 
+ 2c5r + 303X 


famed “1 
where Cor Cyr Cor C3 are unknown coefficients. These equa- 


tions may be expressed as 


c 
w=pPpe 
$= h*c 
where 
p° = Louie, I r°) 
hé —smiseyn ce, len 3r°) 
c = col(c., Cyr Cos C3) P 


20 


If w and $¢? are evaluated at the ends of the element, 


sults may be expressed as 


where 








0 “27 302 


ee 


4 


the re- 


and the subscripts denote evaluation at r=0, and r=%. By 


inverting oe we may solve for c, 


c = Cd 
where 

1 oe ‘i 

lo 16 Oe 

< | -3 -2 22k | 

eo ae 

wan 2 

a2 g2 23g?) 





Bending Energy 


The bending energy, U, of the element is given by 


(13) 


where EI is the flexural rigidity for in-plane bending. 


2 
Now dw 
— ~ = 2c, + 6c,r = t 
dr Z 3 gc 
where g* = owl, OF, “2, Sr) 
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Since EI 1s a scalar constant, the equation may be rewritten 


since c and en are constants, they may be removed from the 
integral. Substituting equation (12) for c and eo we may 
write 
— Lyte: 
U= 7d K'd ae 
where 


K' = bending stiffness = ERIC’ ((aa*ax)c . (145) 
0 


It can be shown by carrying out the required integration and 


matrix multiplications that 


12 62 = 1% 62 
62 49? - 627 Vo sie 

K' = (EI/2%) (16) 
= 2 -62 ee -62 
62 De - 62 49? 


Kinetic Eneuop 


The kinetic energy, T, of the element is given by 
i Q 
tT = 40 war + see 4 }ar (17) 
2 2A 
0 0 
where 
u = mass per unit length 
A = cross sectional area 
I = second moment of area A about an axis through 
the centerline of the element and perpendicular 
to both the w and r axes, 
The first term is the translation energy, T , and the second 
c 


term is the rotational energy, T Thus 


Q 
nog lt -2 
T, = gu dr 
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and 


Q 
eli a @ 32 
r= 3 9 dar 
0 
Now 
w= @, tome & aor’ + cor® 
6 ik D 3 
= pie - p> ce 
Again noting that C and C° do not depend on ug, 
written 
sii SM et 
Ty, - 7a Md 
where 
. eae . Caialte 
M. = uC \ BE Se 
Using ' 292 
od = C7) + 2C5X + 3c4r 
= gh ells aca, 
by a Similar process 
tT. = satu'a 
ie 2—- —r— 
where 
M! - 4 ct(§s ‘ar)< 
—Y A — - a = 
0 
After evaluation 
; 156 222 54 seeps 
228 4y? 132 =o 
My = (uk/420) 
54 132 156 ~22x 
- 1 - 39% - 222 42? 
and 
i" 36 ao eat 32 | 
32 497 —- 32 - 4? 
Mi = (uI/30AQ) 
— 36 ~ 32 36 —- 32 
32 - og = By 42? 
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Bi, 


2 


may be re- 


(18) 


(19) 


(20) 


(21) 


(22) 


(25 


If equations (18) and (20) are recombined, then 


T = 5a'M'd (24) 
where 
M M. + Me (25) 


and M' is called the consistent mass matrix for bending. 


Axial Displacement 

Up to this point there has been no mention of a dis- 
placement in the axial direction. Since the element has no 
axial flexibility, any displacement in the axial direction 
will be displacement of the element as a rigid body. This 
displacement will be the same at either end of the element. 
The axial displacement is represented by u, where U=U =U, - 
In fact, as will be shown later, the axial displacement of 
an entire straight member made up of several elements may 
be represented by a single symbol. 

The axial displacement should be included in the bending 
and kinetic energy equations. The u has no effect on the 


bending energy or rotational energy. It does contribute to 


the kinetic energy of translation. This contribution is 


= ae. Equations (14) and (24) may be rewritten 
U = = v Ky (26) 
where 
= col(u, mo! Por Mg 5) 
POO Che OF 
Q;! 
K = |0! K' (29) 
0; 
Oe 


(30) 


Corner Transformation 

The combination of the element kinetic and bending en- 
ergy to form the total bending energy and kinetic energy of 
the structure is not difficult, except at a corner, Figure 
2a shows two elements, one to the left and one to the right 
of the corner. A subscripting notation has been added. A 
single subscript is used to denote the end or node at which 
the quantity is evaluated. A double subscript is used to 
denote quantities which are the same at each end of an el- 
ement, and for vectors and matrices which pertain to the 


element. Thus Yap col(u,,, Wye oy: Wor do), and so on. 





Wo 2 
mig 
O, 
“2 
(a) Elements Adjacent to (b) Coordinates at 
Corner Node 2 


Figure 2 


om 


The displacements at node 2 are with reference to two 
different local coordinate systems, and cannot be directly 
combined (Figure 2b). The prime is used to denote the dis- 
placements at node 2 which are in the local coordinate system 
of the element to the right of the corner. There are only 
two independent displacements at node 2. It is thus neces- 
Sary to express the four displacements at this node in terms 
of only two displacements. The axial displacements of the 
two elements are chosen. This choice is desirable since the 
axial displacement of all parts of a straight member is the 
same. Further, when end conditions are applied, the axial 
displacements at the ends are some of the quantities which 
are known. 

Letting 6 = cos a and 7 = Sin a, from figure 2b we 
can write 

Our = oor igeliae 7. 
ee ee 


Using the above equations, Wo and Wo, May be expressed in 


terms of Uy > and U513 as: 
“— = —(B/n) ayo 7+ (L/N) 4513 (32) 


This transformation may be put into a convenient matrix 
Mocactton.  Forsihe clement to the left of the corner, this 
transformation is written 


V = Rv* 


t 7, tot 
Vi2 - ia 


32 


where 


i 
| 


| 
t 
: 
Ee 


Es 


po BEL Oy) i ale >>) 
1 0 0 0 0 

0 i 0 0 0 

0 0 1 0 0 
= 

n n | 
oko feeds 


The kinetic and potential energy of the element may there- 


fore be written 


ie) 

Ui 7 
For the element to the 
tion Ls 

o02 a 

= 
where 


(34) 


(a 
= E 9 0 0 
oe 1. 
o> 1 
~<a“ oo oc 1. 


oS. 


The kinetic and bending energy for the element to the right 


of the corner are written 


1 ok t iy ec tok 
To13 = 7¥513) (R') M5, 3R'v5.3 
(36) 
acl, * iE 1, c terk 
Uor3 = 76513) (R') K513R'V5 13 


The kinetic and bending energies in equations (34) and (36) 


can now be combined. 


Example of Method 





(a) Two Member Frame (b) Frame Divided 
Into Elements 


Figure 3 


To illustrate the method, a frame consisting of two 
straight members intersecting at an acute angle, a, is used 
(Figure 3a). Both ends of the structure are fixed. 

Each member is divided into 2 elements (Figure 3b). 

Each of the two members need not have been divided into the 
same number of elements. The nodes at the ends and between 
elements are numbered for identification. The elemental mass 
and stiffness matrices are determined and the following set 
of equations written: 


U (37) 


_ et : it 
D> = Woe ie = 821081212 
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sitar 2’ 

To3 = DM2sMo3Mo3 + Yo3 = 5Y%3ko3%o3 (38) 
ge : “yr 

T34 = DZ3aiaa%3q Ugg = 9% 3.8 ayNs, (39) 
1-t £ 

Ta5 = ZYagius¥a5 « Yay = GMg5kgs5Mas (40) 


The corner transformation must be performed at node 3. 


Equations (38) and (39) become 


A See wey te 
Tyga gute RM RW 7 Ua eer (Voge eee ee 
T,,,= £(7*,,) "(R')  ,, RIV 
ky ey! © 
(42) 
a Se lank 
Dawa Bava» glee Bed Kets yim 


The elemental kinetic and bending energy equations are com- 


bined and are written 


T= =v Vv 
a 2—a—a—a 
(43) 
U_ = Loty V 
a 2—a—a—a 
where 
ve = col(u, 3, eS Pi: bo P51 be P2, War Par Wes >.) 
M,= assembled mass matrix 
K, = assembled stiffmess matriz 


and the subscript a denotes assembled quantities. The assem- 


bled matrices are formed by placing the elemental matrices 
in their correct positions in the larger framework of the 


assembled matrix, and summing all the overlapping terms.“ 


a 


The Equations of Motion 


Using the Lagrangian equation of motion, 


<a /at\_ 9t , au _ 
aul 3 [a 

OV _ — 
equation (43) may be written 


NY SER (4%) 


where gq is thesgeneralized, forcesvectom., The cllemenibqawohac 
are the generalized forces and moments corresponding to the 
elements of the displacement vector, Vo: This equation is 
general in that no restriction has been placed on the gen- 
eralized force vector. 
End Conditions 

We now restrict further development to free vibration 
only. Since in Figure 3 both ends are fixed and there is no 
external loading, the only non=zeromforces in the” vector @ 
are the forces required to satisfy end constraints. 

For both ends fixed, the following is known about the 


structure: 


Tepe | ee a PO 


The order of the matrices in equation (44) is now reduced by 
an amount equal to the number of constraint conditions, in 
this case six. The rows and columns corresponding to these 
conditions are deleted and the equation rewritten 


M"y" + K"v" = 0 (45) 
ao el 
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where 


Ae = col(w., Po O37 War >.) 
KA = a 5x5 reduced matrix 
Ms = a 5x5 reduced matrix. 


Eigenvalue so ieien 


Assuming the vibrational motion is sinusoidal, the 
equation becomes 


K"y" = weM"y" (46) 
— eet —A—a 


where w is the natural circular frequency. This equation is 
now ina form suitable for eigenvalue solution. The solution 
1S carried out in the same manner as described on pages 24 

to 25 of the lumped parameter development. This solution 
will yield the first five natural frequencies (since matrices 
are 5x5), the transverse displacements at the 2nd and 4th 
node, and the slopes at nodes 2, 3, 4. The displacements are 
known at the ends from the end conditions, and through use of 
equation (31) or (32) the displacement at node 3 may be ob- 
tained. These displacements are all referred to local coor- 
dinate systems. They are transformed to a Cartesian (x,y,zZ) 


system through the use of the following equations: 


= —w. in ite ol COS ‘ 
Bs Ww. Sl Va ; Ys 
= . GeoseVect+ we sSaim.y. 
Me Ww. S Yq ; Vs 
where 
y. = angle between the x axis and axial 


direction of member i 


oY 


iioplacement, in x direction 
displacement in y direction 


node number, j= 


member number, 


ac 


1,2,3,4,5 


i=1,2 
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Figure 4 


Gable Bent Used in Analyses 


Figure 4 shows the structure used for comparison of the 
models. The members are all steel wide-flanged I-beams of 
the type shown on the figure. Data for these I-beams is 
taken from reference 5. 

Transfer Matrix Method, 

The transfer matrix method is widely used for the deter- 
Mination of the principal free vibrational modes of framed 
structures and will be used as the standard for comparing 
the finite-element and lumped-parameter-models. The devel- 
opment and use of the transfer matrix method are given de- 
tailed treatment in reference 1. A computer program to an- 


alyze similar structures was developed in 1964 by Fink in 


ow 


his master's thesis at the Naval Postgraduate School. ° The 
“exact values" used to compare the models are generated 
through the use of a transfer matrix computer program devel- 
oped by the author. This program is included as Appendix D. 
The first five undamped natural frequencies and corre- 
sponding mode shapes are calculated using this program for 
four end conditions: clamped-clamped, clamped-simply sup- 
ported, pinned-pinned, and pinned-simply supported. The 
frequencies are given in Table I and the mode shapes are 
listed in Appendix A, pages 50, 51, 52, 53. In Figure 5, 
the mode shapes for the clamped-clamped end condition are 
plotted. The displacements are greatly exaggerated so that 
the mode shape can be more readily visualized. Only the 
first five frequencies are tabulated, since after this the 
computer round-off error begins to affect the mode shape. 
The frequencies have become large enough that, even though 
the convergence routine has estimated them to nine decimal 
places, the frequency determinant is not zero. Thus when 
the state vector of the left end is passed through the trans- 
fer matrices, the accumulation of round-off error is evident 
at the right end. The beginnings of this error can be seen 
in the computer output for the clamped-clamped case (page 50). 
In the fifth mode, the x displacement is 0.0014 at the right 
end, but should be zero since the end is fixed. This sus- 
ceptibility of the transfer matrix method to round-off error 
is a major reason why other computerized models are sought. 
The structure was turned around and the frequencies 
and mode shapes calculated for the pinned-pinned and 
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clamped-clamped end conditions. The frequencies are un- 
changed, but the fifth mode shape differs slightly in the 
fourth decimal place. This difference is due to the accu- 
mulated round-off error. The outputs for these runs are 
listed in Appendix D on page 112. 

Lumped-Parameter-Model 

An analysis of the same gable bent using the lumped- 
parameter-model is next carried out. Since both ends clamped 
is the most severely constrained condition, this end condi- 
tion will cause the structure to have the fewest degrees of 
freedom. Using this end condition, it is expected that cal- 
culation of the natural frequencies and displacements will 
require the greatest number of elements. Therefore, by 
successive trials, sufficient subdividing is carried out 
such that the error in the primary frequency is less than 13% 
and the error in the 5th frequency is less than 3%. To meet 
these criteria requires 38 subdivisions. Using the same num- 
ber of subdivisions, the rest of the end conditions are run 
and the results given in Table I. The pinned-pinned end con- 
dition is seen to lead to a slightly greater per cent dif- 
ference than the clamped-clamped condition. 

A comparison of the mode shapes for this model with 
those of the standard is made by comparing the outputs on 
pages 50, 51, 52, and 53 in Appendix A. The mode shapes agree 
to two decimal places for the first four modes, and to only 
one decimal place in the fifth mode. 

The error due to the lumping of the parameters decreases 


Slowly as the number of subdivisions increases. By comparing 
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Plot of Mode Shapes for a Gable 


Pent With Both Ends Fixed. 


lst mode 


2nd mode 
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3rd mode 





Figure 5 
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page 76 to the values in Table I, the small improvement in 
modal frequencies resulting from increasing the number of sub- 
divisions can be seen. There is a similar improvement in mode 
shapes. If an accurate estimate of frequency and displace- 
ment were desired, a very high degree of discretization would 
be needed. The computer time necessary to solve for these 
values would be large. This is due to the eigenvalue routine 
where the time for solution varies as the cube of the number 
of subdivisions. 

The frame was turned around and values obtained for the 
pinned-pinned and clamped-clamped end conditions. These fre- 
quencies and displacements (page 75) are identical to those 
calculated in Table I. This duplication of results adds a 
measure of confidence to the integrity of the solution using 
this model. 

Finite Element Model 

The frame is again analyzed, this time using the finite- 
element-model. The analysis is started using both ends 
clamped. To meet the criteria of less than 1% difference 
in the primary frequency and less than 3% difference in the 
5th frequency requires 12 elements. The frequencies for all 
end conditions using 12 elements are listed in Table I. 

The mode shapes obtained by this method are in excellent 
agreement with those of the standard. A comparison of the 
mode shapes from this method to those of the standard (Appen- 
dix A, pages 50, 51, 52, and 53) shows agreement to three 


decimal places. 
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TABLE I 


Summary of Modal Frequencies 


Clamped-Clamped End Condition 


Transfer Matrix | Lumped-Parameter Finite-Element 
Sections:6-13- 13- 6 | Sections: ce 4- os Z 


Frequenc 
25 DA 


3 599 
79.647 
150.444 
2.50... 


Transfer Matrix | Lumped-Parameter Finite-Element 
Sections :6-13- 13- 6 | Sections:2-4- a Z 





Pinned-Pinned End Condition 


ode Transfer Matrix | Lumped-Parameter Finite-Element 
umber Sections :6-13-13-6 Sees aa 





Frequenc Frequency |% dif. 
I 15.257 15.285 | +0.182 15.256 =. 008 
2 25.365 25.365 0.00 25.364 |-0.006 
3 75.652 75.162 | -0.65 75.667, |+0. 0K 
4 149.662 M7 Syed 39 149.820 |+0.110 
| 5 277.069 | 269.188 | -2.84 278.160 1+0.390 


Pinned-Simple End Condition 


ode Transfer Matrix | Lumped-Parameter Finite-Element 
Number Sections: - a 13- 6 | Sections: 2-4- 4- 2g 
Frequenc 
9.742 





Ae 


As the number of elements used increases, the differences 
between the calculated values and the standard show a marked 
decrease. The increasing of the subdivisions from 8 to 12 
decreases the per cent difference by a factor of 10. The 
results for the 8 element structure are on page 96 of Appen- 
dix C. This marked decrease in the per cent difference il- 
lustrates the fact that this model will give good results 
with little subdividing. 

It is known that the finite-element-method is equivalent 
to the Rayleigh-Ritz procedure for finding natural frequen- 
cies.> Accordingly the estimate for the frequencies should 
represent an upper bound. The values in Table I do not 
agree with this. In all cases, the first and second frequen- 
cies from the finite-element-model are slightly lower than 
those of the transfer matrix method. An attempt to find a 
reason for this discrepancy in the programming of the models 
was unsuccessful. 

The pinned-pinned and clamped-clamped end conditions 
were rerun with the structure reversed. The frequencies and 
mode shapes were identical to those previously calculated 
(Appendix C, page 95). 

Using the finite-element-model with 12 elements, fre- 
quencies and mode shapes were calculated with rotatory iner- 
tia included. The results are given in Table II. These 
results show that rotatory inertia has little effect on the 
frequencies. The maximum difference between frequencies 


with and without rotatory inertia is 0.3%. 
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Effect of Rotatory Inertia on Frequency 


End 


TAB LE 1 


| Frequency (rad/sec) 

Condition meted: 3 5 

Clamped | W/O’ (25.752 | 34.679 | 79.647 , 150.444 280.597 
| “With 25.748 | 34.675 | 79.585 | 150.207 | 279.725 
‘Clamped | % dif. | 0.015 | 0.012} 0.078] 0.158! 0.312 
Clamped W/O 14.764 | 34.236 | 75.065 | 123.283 | 152.752 
: With 14 34.232 | 75.008 | 123.147 152.516 
Simple 63 dif. | 0.013 | 0.011 | 0.076 0.110 i 0.155 








| 
| 
762 : 
{ 
| 








| 
Pinned | W/O 15.256 25.364 | 75.667 149.820 (278.160 | 
| With 15.255 {25.359 | 75.607 149.583 277.296 | 
Pinned | % dif. _ 0.006 | 0.020! 0.078; 0.158; 0.310 
Pinned | W/O | 9.741 | 24.417 | 72.502 | 122.956 (152.412 | 
With | 9.740 , 24.413 | 72.415 | 122.819 | 152.177 
Simple | 3 dif. | 0.001} 0.016 | 0.079, 0.112} 0.159 | 


Comparison of Models 


Both the lumped-parameter and finite-element-models 
lend themselves favorably to computer analysis. The param- 
eters needed for analysis are easily derived and fit conven- 
iently into a matrix form which is easily programmed. Each 
model has one feature which may cause some difficulty. These 
are the application of end conditions in the lumped- 
parameter-model and the corner transformation in the finite- 
element-model. Once understood, however, both of these can 
be programmed. 

The inertial and stiffness matrices resulting from these 
models are quite different. The matrices from the finite- 
element-model are sparse. The matrices from the lumped- 
parameter-model do not have this property. The eigenvalue 
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solution will therefore be easier and less time consuming for 
the finite-element-model since the matrices are already par- 
tially diagonalized. 

The order of the resulting matrices is also different. 
The order of the finite-element-matrices is 2N+3 where N is 
the number of elements. The order of the lumped-parameter- 
matrices is N+l. The finite-element-model using 12 elements 
or 27 degrees of freedom yields much better solutions than 
the lumped-parameter-model with 38 segments or 39 degrees of 
freedom. This significant difference in the order of the 
matrices will save computation time also. 
Conclusions 
1. Both models are conveniently treated by digital computer 
methods and could be used for the analysis of frames. 
2. Both models can include rotatory inertia. 
3. Both models provide a basis for finding dynamic response 
BE the structure to arbitrary time-varying external loads. 
4. A comparison of the results of an analysis of a gable 
bent using the models shows that the finite-element-model 
gives much better results. It requires 12 fewer degrees of 
freedom and its average difference from the standard is one- 
tenth that of the lumped-parameter-model. 
5. The transfer matrix method still gives good results for 
frequencies higher than the fifth, but not mode shapes. The 
finite-element-model gives good agreement at these higher 
frequencies, and its mode shapes should -be caccurate <sinte ‘the 
particular round-off error that affects the transfer matrix 


method does not pertain to the finite-element-model. 
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APPENDIX A 











TPANSFER MATPIX METHOO 


x COQRB LIN) Y COOROLINY 


30000000.00000 
0-0 60.00000 
20000000.00000 
0.0 120.00000 
30000000.00000 
180.00000 156.00000 
30000000.00000 
360.00000 192.00000 
30000000.00000 
540.00000 156.00000 
30000000.00000 
720.00000 120.00000 
30000000.00000 
720.00000 60.00000 
30000000 .00000 
720.00000 0.0 
CLAMPEO CLAMPEO ENO CONOTTION 
ROTATOPY INEPTIA [S NOT I[NCLUOEO 
NATURAL FREQUENCTES€RAO/SEC) 
l ? 3 4 
1 25.754 34.680 79.630 150.284 
NORMALTZEO X OLSPLACEMENTS 
I 2 3 4 = 6 
1 0.0 ~0.0365 ~0.1476 -0.2313 -0.2578 -0.301 
2 0.0 0.1980 0.4195 0.3988 0. 3440 2 346 
a 0.0 0.1604 0.1948 0.0614 0.1813 2042 
5 0.0 0.6216 0.0224 -0.0149 0.0352 2043 
NORMALIZED ¥ OITSPLACEMENTS 
1 2 3 4 5 6 
i 0.0 0.0000 0.0000 0.4183 0.5509 5 3)3h) 
2 0.0 -0.0000 -0.0000 0.1038 0.3777 2391 
3 0.0 -0.0000 -0.0000 0.6670 0.0674 -0.627 
% 0.0 -0.0000 0.0000 0.5262 -0.6017 o4TT 
5 0.0 -0.0000 -0.0000 0.1875 -0.0633 =0.021 


LUMPFOQO PAPAMETER METHNO 


x COpReCING vf CORE Re UN! 


0.0 120.00000 
360.00000 192.00000 
720.00000 120.00000 
720.00000 0.0 


CLAMPED CLAMPEO ENO CONOITION 
FOTATOPY INERTIA TS NOT INCLUDED 


NATUPAL FPEQUENCTES(RAO/ SEC) 


l 2 3 
1 25.921 34.929 79.4456 
NOPMALTZEO KX OTSPLACEMENTS 
1 2 3 4 
1 0.0 0.0352 0.1452 0.2294 0 
Fis 0.0 0.1969 0.4192 0.4000 0 
3 0.0 -0.1613 -0.1984 -0.0653 =(0) 
4 0.0 0.0600 -0.1054 -0.2119 0 
5 0.0 0.6208 0.0247 -0.0160 0 
NORMALTZEO ¥Y OTSPLACEMENTS 
1 2 3 4 
l 0.0 0.0 0.0 -0.4202 S10, 
2 0.0 0.0 0.0 0.0954 0 
3 0.0 0.0 0.0 -0.6654 =) 
% 0.0 0.0 0.0 0.5324 -0 
5 JO ORO 0.0 0.2034 -0 


FINITE FLEMFNT METHOO 


X CAORDOCIN MCR 
Pia ) / (¢ peg iN 


0.0 120.09000 
360.09000 192.00300 
721.00000 120.00000 
720.00000 0.0 


CLAMPED CLAMPEO ENO CONDITION 
KIJTATORY INERTIA [TS NOT INCLUDED 


NATURAL FREQUENCLESCRAO/SSFC) 


1 ? 3 
1 BD aU Gi¢2 34.679 79.4547 
NIKMALIZEY X NTSPLACFMENTS 
1 2 4 4 
i 0.0 -0.0365 -0.1474 =o erie = 
2 0.9 90.1580 0.4196 0.3688 0 
3 0.0 -O0.1494% -0.1944 -0.0612 (0, 
4 0.0 -J.0587 0.1077 9.2130 10 
5 0.0 Sl ging ihs' -0.02272 0.0153 ~0 
NORMAL ET/JFO Y OTSPLACEMENTS 
1 2 3 4 
1 9.0 1.9 0.9 0.4183 0 
? 0.U 0.0 0.0 0.1034 0 
3 0.0 0.0 0.0 =-0.45467)] -0 
4 0.0 0.0 0.0 -0.5263 0 
= 0.0 0.0 0.0 ~0.1872 0 


4.00000 


13.00000 
13.00000 


6.00000 


4 
148.26? 


2.00000 
4.00000 
4.00000 
2.00000 


MOOULUS OF FLASTICITY(PST) 


NUMBEP OF SECTIANS 


6 
58 0.2999 
44 0. 3483 
45 -0.0455 
34 0.2784 
79 0.0464 


NUMBEP OF SECTIONS 


30 


0.28300 
0.78300 
0.28300 
0.28300 
0.28309 
0.78300 
0.28300 
0.29200 


5 
279.474 


0.28300 
0.28309 
0.28300 
9.29300 


5 
271.688 


0.28300 
0.22300 
0.278300 
9.29300 


SPECIFIC wWrtlL8/Ctl 


SPECIFIC WT(L8B/cU 


SPFCLFIC WT(L&/CU IN) 


TN} 


IN) ARFACSO ITN) 


5.8RNN0 
5 .ARONN 
29.40000 
29.40000 
29.40000 
29.40000 
4.62000 
4.67009 


AREATSO IN) 
5.28000 
29.40000 
29.409000 
4.42000 


APFA€SO IM) 
5.88000 
29.40000 
29.49000 
4.62000 


MOMENT( TA 4) 
A2 50900 
27.50000 

625.90N00 
625. 99000 
é75.09900 
A25.90000 
30.3000900 
30.20000 


MOMENTC IM 4) 
82.50000 
6?5.90000 
625.00090 
329.30000 


MOMENTCIN 4) 
&2.500900 
£75. 97990G 
625.90900 
20.20000 


TRANSFER MATRIX METHOO 
x CAQRACING 


0.0 
0.0 


140 .00000 
360 .00000 
540.00000 
720.00000 
720.00000 
72000000 
CLAMPEO SIMPLE ENO CONDITION 

ROTATORY INERTIA IS NOT INCLUOED 


Y COMROLIN) 
6° 


NATURAL FREQUENCIES(RAO/SEC) 


1 


NORMALIZEO X OLSPEACEMENTS 


NEW 


SI ys NS me 


l 
14,765 


Z 
0.0000 
-0.0000 
-0.0000 
-0.0000 
-0.0000 


LUMPEO PARAMETER METHOO 
x are aN 


MODULUS OF ELASTICITY¢(PSI) 


429000000 .00000 
460.00000 
30000000.00000 
120.00000 
30900000.00000 
154.00000 
29000000.00000 
192.00000 
30000000 .00000 
156.00000 
19000900.00000 
120.00000 
30000000.00000 
40.00000 
39000000.00900 
0.0 
2 3 4 
34.228 USO 52 123.734 
~ 5 § 
-0.1537 -0.70845 -0.2302 -0.2538 
0.40469 0.37481 0.3217 0.3200 
0.0960 0.0184 0.0800 -0.0027 
—mUe -0. 1031 -0.0016 0. 1011 
3 ~ 5 4 
0.0000 0.2754 0.358273 0.244% 
-0.0000 0.0463 -0.03457 -0.0171 
0.0000 0.2339 -0.2734 0.2403 


Y CAMROCIN) 
0.0 


NUMBER NF SECTIONS 


4.00000 
0.0 120.900000 
13.00000 
360.00000 192.00000 
13.00000 
720.00000 120.00000 
6.00000 
720.00000 0.0 
CLAMPEO SIMPLE ENO CONOITION 
ROTATORY INERTIA [S NOT INCLUOED 
NATURAL FREQUENCIES(RAO/SEC) 
1 2 3 4 
1 14.769 34.482 74.760 123.969 
NORMALI7ZEO X NISPLACEMENTS 
1 2 3 5 4 
1 0.0 0.0500 0.1526 0.2 0.2293 0.2529 
Zz 0.0 0.1943 0.4067 0.3793 0.3222 0.3214 
3 0.0 -0.0954 -0.0992 -0.02190 -0.0827 0.0019 
5 0.0 -0.0213 0.0507 2093 0.9022 -0.0904 
NCRMALIZEO Y OISPLACEMENTS 
1 2 3 & 5 4 
1 0.0 0.0 0.0 -0.2T4AT -0. 3535 -0.2659 
2 0.0 0.0 0.0 0.1371 0.4225 0.4195 
3 0.90 0.0 0.0 -0.3905 -0.0823 0.3363 
4 0.0 0.U 0.0 0.0483 -0.0355 -0.0140 
5 0.0 0.0 0.0 -“0.2117 0.2423 -0. 2208 


PPSITE ELEMENT “ETHNO 
xX CIOROCIN) 
0.0 


Y COOROCIN) 
0.0 


NUMBER OF SECTINNS 


?.0000 
0.0 120.00000 e 
4. 
359.00000 192.09000 ae 
4.00009 
720.00000 120.900900 om 
> 
720.000U0 0.9 ee 
CLAMFEO SIMPLE ENO CONDITION 
ROTATNRY INERTIA IS NOT INCLUDEO 
NATURAL FIEQUENCIFS(RAO/SFC) 
l 2 4 4 
l 142 746 344234 75.065 123,783 
NAKMALIZED x DISPLACEMENTS 
1 ? 3 4 5 6 
1 0.0 -0.0504 -0.1537 -0.208R -0.2302 -0.2538 
2 0.0 0.1953 0. 4064 0.3781 0.3217 0.3200 
3 0.9 -0.0839 -0:0959 -0:.U01H4  -0.0800 939027 
4 9.0 “1.9108 -0.0082 0.0011  -0.0153  =9.0190 
5 9.0 0.0234 -0.0564 -021033 -0.0016 0.1913 
NOKMALIZED Y DISPLACEMENTS 
1 > 3 4 5 2 
1 9.0 0.0 0.9 0.2754 0.3973 0.2444 
2 0.0 0.0 v.0 0.1438 0.4256 On 4 172 
3 0.0 9.9 0.0 -0.3845 -9.0797 0.33237 
4 0.0 0.9 0.0 -0.0463 9.0357 9.0171 
5 0.9 9.0 0.0 0.2345 -0.2739 0.2407 


SPECIFIC WT(LS/CU 


0.78300 
0.28300 
0.28300 
0.28300 
9.25300 
0.28300 
0.242300 
0.29300 


0.29309 
0.28300 
0.293090 
0.25300 


0.25300 
9%. 28300 
0. 278300 
0.278700 


5 
152.752 


-0.0000 
0.00900 


-0.0000 


IN) 


SPECIFIC wWT(LB/CU IN) 


A 


2 
? 


re 
A 


RFA( SQ TN) 


§.A8NN0 
5.88000 
9,40009 
7.490000 
9.40NNN 
9,40000 
4.£70090 
4.47000 


AQFACSO TN) 


5.98000 


29,40000 
29.40N00 


4.6700N 


0.09090 


-0.0090 
-9.0000 
-0.0000 
-0.0000 


SP9FCIFIC WTHt(LASCIU IN) 


AQFAISO IN) 
5.AR000 


?9,4N000 


29.40000 


4 .A2N0N 


MAMENTE TN 6) 
a>.snnan 

R72 ,F NAAN 
£28 VANAN 
425.99009N 
425,00N99 
APRLNIANN 
aN. 390NnNn 


30.79009% 


MOMENTO AD 4) 
G7 SAIAA 
675.99099 
625.390090 
49.3099 


MOMENTOITN 4) 
2>.8qnnn 
A25.9900NN 
426.7900 
39.70999 


TRANSFER MATRIX METHOD 
x CAQROCINI 


30000000.00000 
0.0 60.00000 
30000000.00000 
OO 120.00000 
30000000.00000 
180.00000 156.00000 
30000000.00000 
360.00000 192 .00000 
30000000.00000 
540.00000 156.00000 
40000000.00000 
720 .-00000 120.00000 
420000000.000090 
720.00000 60.00000 
30000000.00000 
720.00000 0.0 
PINNED PINNED END CONOITION 
RDTATGRY INERTIA IS NOT INCLUNED 
NATURAL FREQUENCIES(RAD/SEC) 
1 2 3 4 
1 15.257 25.365 75.652 149.662 
NORMALIZED X DISPLACEMENTS 
1 2 3 4 5 6 
1 0.0 0.2170 0.3622 0.3957 0. 3960 0.4152 
2 0-0 0.2132 0.2569 Oreletar 0.1257 0.0927 
3 0.0 0.1993 0.1320 -0.0020 0.1233 -0.0153 
4 0.0 0.0961 -0.1133 -0.2159 0.0078 0.2238 
5 0.0 0.6120 -0.0077 -0.0229 -0.0013 0.0061 
NORMALIZED Y OISPLACEMENTS 
1 2 3 4 5 6 
1 0.0 -0.0000 -0.0000 -0.1678 -0.1693 -0.0733 
2 0.0 ~0.0000 -0.0000 0.3961 0.6559 0.4906 
3 0.0 -0.0000 -0.0000 0.6702 0.0438 -0.6493 
o 0.0 -0-0000 0.0000 0.5130 -0.6055 0.4743 
5) 0.0 -0.0000 0.0000 0.0758 -0.0323 0.0044 


LUMPED PARAMETER METHOD 
x COBRDCIND 


0.0 
360.00000 
720.00000 
720.00000 


PINNEO PINNED END CONOITION 
IS NOT INCLUOEO 


ROTATORY INERTIA 


Y Oe OLN 
~ O. 


Y COGRRTINI 


120.00000 

192.00000 

120.00000 
0.0 


NATURAL FREQUENCTES(RAO/SEC}) 


l 
I 15.285 


NORMALIZED X DISPLACEMENTS 


1 2 
1 0.0 -0.2149 
2 0.0 -0.2140 
3 0.0 -0.2003 
4 0.0 -0.0991 
5 0.0 -0.6146 


NORMALIZEO ¥ DISPLACFMENTS 


1 2 
1 0.0 0.0 
é 0.0 0.0 
3 0.0 0.0 
4 0.0 0.9 
5 0.0 0.0 


FINITE cLEMENT METHNO 
xX COQRRCINI 


0.0 
360.00000 
720.00000 
720.00000 


PINNED PINNED ENO CONDITITCN 
1S NOT INCLUDED 


RMTATCRY INERTIA 


2 3 

25.365 75.1462 

3 % 
-0. 3618 -0.3956 -0 
~0.2574 -0.1783 -0 
o)5 haves 0.0013 10, 
0.1115 0.2152 -0 
0.0081 0.0243 0 

a 4 
0.0 0.1689 0 
0.0 -0.3953 -0 
0.0 -0.6698 (0, 
0.0 ~0.5178 0 
0.0 -0.0807 0 


Y COOROIIND 


NATURAL FREQUENCIESIRAO/SEC) 


1 
1 15.256 


NORKMALIZEO X OISPLACEMFENTS 


1 2 
1 0.9 U.-2170 
2 1.0 0.2132 
3 0.0 0.1993 
& 0.0 -0.0941 
5 0.uU -0.6119 


NOKRMALIZEO Y OTSOPLACEMENTS 
] 2 


SH wrere 
oOcecoo 
@eeese 
SCooce 


MDDULUS DF ELASTICITY¢(PSI) 


NUMRER DF SECTISNS 


6.00000 
13.00000 
13.00000 

6.00000 

4 

5 6 
3958 -0.4151 
1264 -0.0936 
1238 0.0148 
0081 -0.2234 
0015 -0.0062 

3) a 
1698 0.0733 
6550 “0.4913 
0444 0.64Rf 8 
5988 “0.4775 
0332 -0.0054 


NUMBER OF SECTIONS 


2.00000 
120.00000 
4.00000 
192.00000 
4.000900 
120.00000 
2.00000 
0.0 
2 3 4 
25.7264 75.647 149.820 
3 4 5 6 
0.3622 0.3957 9.3960 0.4152 
0.7569 O.1777 0.1257 0.0927 
0.1319 -0.9021 0.1232 -0.0154 
0.0078 0.0229 0.0013 -0.0960 
3 4 5 6 
0.0 ~0.1678 -0.1693 -0.0733 
0.0 0.3961 0.6559 0.4906 
0.0 -0.5130 0.6054 “0.4744 
0.0 -0.0754 On osz2 -0.0044 


a2 


0.28300 
0.28300 
0. 28300 
0.28300 
0.28300 
0.28300 
0.28300 
0.28300 


5 
277.04 


0.28300 
0.28300 
0.2°200 
J.28300 


S 
269.18 


-0.0000 = 
0.0000 
-0.9000 c 
0.9000 
-~0.9000 < 


0.28300 
0.78309 
0.28309 
0.28300 


5 
278.14 


SPFCIEIC WTtLA/CU IN) 


9 


SPECIFIC WT{LA/CU IN} 


8 

a 
0.2987 
0.0489 
0.2048 
0.1024 
0.7845 


SPECTFIC WY(LA/CI IN) 


0 


ARFEA(SQ IN) 
5.88009 
5.88000 

79.40000 

29.40000 

29.40000 

29.40000 
4.62000 
4.62000 


ARFAISQO 1M) 
5.88000 
29.40000 
29.40000 
4.62000 


-0. oon0 


AREA( SOQ IN) 
5.88000 
29.40000 
29,49000 
4.42000 


MOMENTE TN 4) 
R2.50000 
8?,500N0 

625,00000 
625.0000) 
675.90000 
625.00000 
30.30000 
30.30000 


MOMENTOCIN 4) 
R2.£00N0 
425.00000 
4625.00090 
20,20900 


MIMENTO(IN 4) 
8?.50000 
425.30000 
6?5,.90000 
30.320000 


TRANSFER MATRIX METHOD 


x COBRRCINY V/ COOROLING MODULUS OF ELASTICITYEPSI} SPECIFIC WTILA/TU IN} AREATSO IN) MOMENTOCIN 4} 
° : 30000000.00000 0.28309 5.AROOO R?2.50N90 
0.0 60.00000 
30000000.00000 0. 2R300 &.AROOO A>. 5NONO 
0.0 120.00000 
390000000.00000 0.28300 29.40000 625.9NNAD 
180 .00000 156.00000 
30000000.00000 0.78300 29.40000 675.90000 
360.00000 192.00000 
30000000.00000 0.28300 29.40000 4256.900N9 
540.00000 156.00000 
30000000.00000 0. 2R300 29.40000 475.9nNNnON 
720.00000 120.00000 
30000000.00000 0.28300 4242000 Siiies SA a ade 
720 .00000 60.00000 
30000000.00000 0.28300 4e.67NNN 3N.3000Nn 
720 .00000 0.0 
PINNED SIMPLE END CONOITION 
ROTATORY INERTLA IS NOT LNCLUOED 
NATURAL FREQUENCTES(RAO/ SEC) 
1 2 3 4 5 
1 9.742 24.419 72.489 122.907 152.252 
NORMALIZED X OLSPLACEMENTS 
1 2 3 “6 5 6 7 a 9 
1 0.0 -~0.1556 -~0.2682 -0.3028 -0. 3120 -0.3270 0.355R ~0.4092 -7.463R 
Z 0.0 0.2354 0.3043 0.2398 0.1903 0.1668 0.0743 -0.1014 -0.2R30 
4 0.0 0.0152 0.0061 -0.0027 0.0137 0.0183 0.0213 -0.3189 -0.9452 
5 0.0 0.0377 -0.0589 -0.1035 -0.0051 0.0961 0.04ART7 On 257) 0.8532 
NORMALIZED Y DISPLACEMENTS 
1 2 3 4 5 6 7 a a 
1 0.0 0.0000 0.0000 Onle2 0.2190 0.1439 0.0000 0.0090 0.N000 
2 0.0 -0.0000 -0.0000 0.3227 0.5701 024527 -0.0000 -~0.0000 -0.0000 
3 0.0 -0.0000 -0.0000 0.4026 0.N578 -0.3690 0.0000 0.0000 0.N000 
4 0.0 -0.0000 -0.0000 0.044% -~0.0379 -0.0151 0.0000 0.0000 0.0000 
5 0.0 -0.0000 0.0000 0.2230 -0-2689 0.2375 0.90000 9.0000 0.0000 





LUMPED PARAMETER METHOD 


x CORBRETN MOMENTEIN 6} 


Y COQRBIIN) NUMRER SF SECTIONS SPECIFIC WIt{LA/Cty IND AQFALSO IN) 


6.00000 0.2R300 5.88000 a?.S000N 
° 120.00000 
one 13.00000 0.28300 29.40000 62?5.nN00NN 
60.00000 192.00000 
5 13.00000 0.28300 29.40000 625.09N00 
720. 120.00000 
ps dag 6.00000 0.28300 4.62000 30. 3400N0 
720.00000 0.0 
PINNED SIMPLE ENO CONDITION 
RDTATORY INERTIA IS NOT INCLUDED 
NATURAL FREQUENCIESIRAD/SEC} 
1 2 3 4 5 
1 9.747 24-419 72.041 123.511 150.555 
NORMALIZED X DISPLACEMENTS 
1 7 a 4 5 os) y a 9 
1 0.0 0.1555 0.2679 0.3026 0.3118 0.3768 0.3557 0.4092 0.4639 
2 0.0 -0.2361 -0.3049 -0.2406 -0.1911 -0.1680 -0.0773 0.1007 0.2R22 
3 0.0 0.1140 0.0714 -0.0101 0.0597 -0.0268 0.0479 pares 0.7370 
4 0.0 0.0156 0.0054 -0.0040 0.0135 0.0193 0.0217 -O0.31R4 -0.9451 
5 0.0 -0.0344 0.053! 0.0933 0.0056 -0.0355 -0.0419 -0.2659 -O.R7T3A7 
NORMALIZED Y DISPLACEMENTS 
1 2 3 4 5 6 i a 9 
1 0.0 0.0 0.0 -0.1740 -N.2195 -0.1444 0.0000 0.0000 0.0000 
2 0.0 0.0 0.0 -0.3216 -0.5691 -0.4534 0.9000 9.0007 0.9000 
3 0.0 0.0 0.0 0.4073 0.0587 -0.3734 0.0000 0.0000 0.0000 
4 0.0 0.0 0.0 0.0470 -0.0409 -O.O1I1R 0.0000 -0.9000 -0.0000 
5 0.0 0.0 0.0 -0.2006 0.2376 -0.2178 0.0005 0.0000 -0.000N 
me SS eee EEE 
FINITE FLEMENT METHOO 
x ee os aa Y CADE RAIUN? NUMRER NE SECTIONS SPECTFIC WHt{l8/Ci) IN} AQFACSO TN} MOMENTE TS 4) 
5 2.00000 0.28300 5.2A000 R?.609NN 
0.0 120.00000 
4.90000 0.28300 29,4000N #78 .009N0 
360.00000 192.00900 
4.00000 0.2R300 29.40009 ¢76.00N10 
720.00000 120.0000 
2.900000 0.7R300 4.267000 30,.30000 
720.00000 0.0 
PINNEO SI{“PLE END CONOTTION 
RCTATGRY INERTIA IES NOT INCLUNEO 
NATURAL FREQUENCIES(IRAD/SEC) 
1 2 a 4 5 
l 9.741 24.417 72.502 122.956 152.412 
NOIWMALTZEN xX OTSPLACEMFENTS 
1 2 u & 5 4 1 a 9 
1 0.0 0.1556 0.2687 0.302R 0.3120 0.3270 0.3558 0.4N92 0.4638 
2 0.0 0.2354 0.3043 0.2398 0.1993 O.144R 0.0TA3 -O.1NIT -0.2A%3%0 
3 0.0 -O.-1119 -0.0699 0.0106 -0.0584 0.0270 ~0.0468 -0.3570 -0.7421 
4 0.0 -0.0152 -0.0061 0.0077 = Om Olea -O.-O1R2 -0.02713 0.3}A49 0.94527 
5 0.9 -0.03 78 0.0590 0.1037 0.0051 -0.0964 -0.04A8 -0.2574 -O.R528 
NOKMALIZED Y DLSPLACEMENTS 
I ? 3 4 5 4 7 q Q 
1 0.0 0.0 0.0 ~0.1732 -0.2190 ~-0.1439 0.0000 0.0 0.0 
2 0.0 0.0 0.0 Owere 0.5702 0.4527 0.00N0 0.0 0.0 
3 0.0 0.0 0.0 -0.4025 -0.0578 0.34690 -0.N000 0.0 0.0 
4 0.0 0.0 0.0 -0.0444 0.0378 0.0152 -0.NNOO 0.0 0.0 
5 0.0 0.n 0.0 -0.2235 0.7694 -0.2379 -0.9009 0.0 0.0 


oS 


FINITE ELEMENT METHOO 


4) 
R?P.509N4 


MOMENT ( FA! 


IN) ARFA(SO IN) 
5.8000 


0.28300 
0.28300 
0.28300 
0.28300 


SPECIFIC WTELB/CU 


NUMBER DOF SECTIONS 
2.00000 
4.00000 
4.00000 
2.00000 


Y COORRCTINI 


Rp Cli 


x COR 


675.79N0N 


29.40000 
29. 40000 
“.62000 


120.00000 
192.00000 
120.00000 


0.0 
360.00000 


625.90099 


30,3279 


72Q.00000 


720.00000 
CLAMPED CLAMPED END CONOITION 


ROTATORY INERTIA 


0-0 


INCLUDED 


iS 


NATURAL FREQUENCIES(RAD/SEC) 


S 
UD CASS 


és 


150.207 


3 
34.675 79.585 


1 
25-748 


NORMALIZEC X OISPLACEMENTS 


ANN TW 


NORMALIZED Y DISPLACEMENTS 


MINION TI 


FINITE ELEMENT METHOO 


MOMENTOCIN 4) 


Y COQRRIIN) NUMBER DF SECTIONS SPECIFIC writ B/C IN) AREA(SO IN) 


x COORBCIND 


82.5059" 
625.9098 
625.00090 


5.PRO0O 
29.40000 
29.40000 


0.28300 
0.28309 
0.28300 
0.28300 


2.00000 
4.00000 
4.00000 


120.00000 


0.0 
360.00000 


72000000 
720.00000 
CLAMFEC SIMPLE ENC CONDITION 


RCTATORY INERTIA IS INCLUDED 


192.00000 
120.00000 


40.3009N 


4.62000 


2.00000 


0.0 


NATURAL FREQUENCTES(RAC/SEC) 


S 
162.516 


4 


123.147 


5 
75.008 


14.762 


NORMALIZED xX OISPLACEMENTS 


HINA SW 


NCRMALIZEC ¥ OILSPLACEMENTS 


ANT SW 
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FINITE ELEMENT METHOO 


a} 


S 

cS 
«Cc 
=- CC 
~— ww 
—- . 
ze 
Ww «& 
> 3 
= 
> 2 
Zz 

Cc 
oc 
Ww © 
—- 
a « 
Ww «+ 
ar w 
<a 
z 
— 
-~ 
Stl 
= 
~ 
c. ©] 
—) (& 
-—- © 
- ~- 
zx « 

A 
Ww e 
- © 
uw 
=) 
Wu 
a 
“" 
VV) 
z 
oe) 
= 
UU Oo 
w & 
wo 

o 
uw © 
CG « 

~ 
a 
iw 
® 
= 
2 
Zz 
CO 
ae 
OO 
= 
Ww 
~ 
z 
co 
ae 
Oo 
SZ) 
1) 
4 


426 .0N000 
62°.30900 
InN, 3NANY) 


29.40000 
29.40000 


0.28309 
0.28300 
0.28300 


4.00000 
4.00000 


120.0Vv0C0 
192.00000 
120.00000 


0.0 
360 -00000 


4.62000 


2.00000 


720.00000 


0.0 


720.00000 
PINNEO PINNEO ENO CONDITION 


ROTATORY INERTIA [LS INCLUOEO 


NATURAL FREQUENCTES(RAO/SEC) 


4 5 
277.296 


149.583 


3 
75.607 


2525259 


1 
15.255 


NORMALIZEC X OTSPLACEMENTS 


NN TUN 


NORMALIZEC Y OITSPLACEMENTS 


mA PU 


S(NTTE ELEMENT METHOO 


MOMENT! 4) 


IN} AREAESO IN) 


SPECIFIC WTI{LA/CY 


NUMBER OF SECTIONS 


Y CORRREIND 


x CONRRCIN) 


R2.59C90 
&7?5.09910 
62£.00000 

30.20000 


5.88000 
29.40000 
29.40900 


0.28300 
0.29200 
0.28300 
0.28200 


2.00000 
4.00000 
4.00000 
2.00000 


120.00000 
192.0000 
12v.00000 


0.0 
360.00000 
720.00000 


4.627000 


0.0 


720.00000 
PINNEO SIMPLE ENO CONDITION 


ROTATCRY INERTIA 


1S INCLUNEO 


NATURAL FREQUENCIES(RAN/SEC) 


152.177 


4 
122.819 


3 
72.445 


24.413 


9. 74) 


NIR“M“ALTZEC X OTSPLACEMENTS 


an FUN 


NCRMALIZEC ¥ CISPLACEMENTS 


me NY OS OS 
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APPENDIX B 


LUMPED-PARAMETER PROGRAM 

General 

Program LUMPED is a double precision IBM FORTRAN IV 
Language’ digital computer program designed to find the 
principal vibrational modes of straight sided plane frames 
of uniform cross section. The program is written for the 
IBM 360/67 facility installed at the Naval Postgraduate 
School. The in-plane undamped natural frequencies, as well 
as normalized displacements, may be found through its use. 
The displacements are calculated for the mid-point and ends 
of each straight member. The analysis is performed by the 
characterization of the structure as a chain of pin-connec- 
ted rigid bars with both mass and flexibility lumped at the 
hinges. The flexibility is represented by torsional springs. 
The model neglects axial and shear deformation, but rotatory 
inertia can be included if the user desires. 

The program 1S dimensioned so that a structure with 
up to nine straight members may be analyzed. These members 
may be divided into a maximum of 79 segments (total). The 
program with this dimensioning requires total storage of 
400K bytes. 

The more severely constrained end of the structure is 
designated the left end. This end must be either fixed or 
pinned. The right hand end is allowed any of the four 


common boundary conditions: fixed, pinned, simply supported, 


56 


or free. The structure may not have any intermediate sup- 
ports or externally applied forces. 
Program Structure 

The program consists of a main body and two subroutines. 
If double precision sine, cosine, arctangent, absolute value, 
and square root routines are not available in the FORTRAN 
library, they must be provided by other means since they are 
required by the program. The main program contains comment 
cards explaining each major section's function, and a large 
leading section giving input format information. 

When the end conditions are applied, the order of the 
matrices is reduced. Rather than rearrange the storage of 
these matrices, the rows and columns corresponding to the 
no-longer-needed equations and rotations are replaced with 
zeroes. In the mass matrix, each diagonal element corre- 
sponding to a deleted row and column is replaced with unity. 
The eigenvalue problem may now be undertaken. The final 
spectral matrix will contain zeroes for the eigenvalues of 
the deleted rows. 

The subroutine JACVAT obtains a solution for the real 
eigenvalue problem of a symmetric matrix. The subroutine 
is based on Jacobi's diagonalization by successive rota- 
pions.” The subroutine will take a matrix up to the maximum 
size of 160 by 160. Care must be taken to insure that the 
entering matrix is symmetric. If any of the off-diagonal 
elements are non-symmetric, the subroutine prints the loca- 


tion of these elements, and then makes them symmetric by 


a7 


averaging them. This output will be printed in the midst 
of the regular output, and is therefore not desirable. 

The subroutine MAPRIN is the output subroutine. This 
subroutine will print a matrix, column by column. It num- 
bers the rows and columns and allows a single line heading 
of 32 spaces. The heading must be in a DATA statement in 
the calling program. The subroutine has six different 
formats which may be used. 

Input 

For program input, the following things are needed: 
1. Coordinates of the ends and the corners of the struc- 
ture. A Cartesian coordinate system is used, with the 
positive x direction to the right and the positive y direc- 
tion up. All coordinates are in inches, and start with 
the more severely constrained (left) end. 

2. Number of segments into which each member is to be 
divided. 
3. Properties of each member. 

a. specific weight, lb/cu.in. 

Did be crs Gen. 

C.. modw@itis=er elasticity, psi. 


d. second moment of the cross sectional area with 
respect to the z axis, in. 


4. End conditions. 
5. Angle at which right hand end condition is applied, 
from +90° to -90°. This angle is of concern only when the 


Eaght jendwismsimplyssupported.  Forwall.othexecases,,! "SUMS 
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taken as zero. The angle is the acute angle between the y 
axis and the line of action of the restraining force. It 

is measured from the y axis and is positive counterclock- 

wise. 

The FORMAT for input data, as well as the make-up of 
the data deck, are listed in the comment section at the 
start of the main program. 

Output 

The program will automatically print for each problem 
and set of end conditions within a problem the following: 
1. Coordinates of the structure and number of segments 
between each set of coordinates. 

2. The properties of each member. These properties are 
the same as those in the input section. 

fe «Ene@econdrtions. 

4. Whether or not rotatory inertia has been included. 

The user has control over the rest of the output. 

The user specifies the number of natural frequencies de- 
sired with card number 220 in the program listing. The 
total number of frequencies requested may not exceed the 
number of degrees of freedom of the model minus the number 
of constraints. The specified number of frequencies will 
be printed in ascending order. The mode shapes correspond- 
ing to the frequencies are printed in two arrays. The u 
components are printed first and the v components second. 


Fach row of these arrays corresponds to a particular 


ys) 


frequency. The frequencies increase going down the columns 
of the arrays. Each mode shape is normalized such that the 


Square root of the sum of the squares is unity. 
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FOR USE OF LUMPED PARAMETER PROGRAM 
NUMBER OF PROBLEMS TO BE WORKED 
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C-----APPLICATION OF END CONDITIONS USING A GAUSS ELIMINATION 
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LUMPEO PARAMETER METHOO 


MIMENTCTN 4) 


ARFACSO IN) 
4.62000 
79.40000 
29.40000 
5.88000 


SPECIFIC WFK{LB/CU IN) 


OF SECTIONS 


6.00000 
13.00000 
13.00000 


NUMBER 


Y COORRCIN) 


Xx COORRIIN) 


30.79000 


62?5.00090 
&25.90000 


0.28300 
0.28300 
0.28300 
0.28300 


120.00000 
192.00000 
120.00000 


0.0 
360 .-00000 
72000000 
7720-00000 
CLAMPEO CLAMPEO ENO CONOITION 


R2.59000 


6.00000 


ROTATORY INERTIA IS NOF INCLUOEO 


NATURAL FREQUENCTES(RAO/SEC) 


5 
271.688 


3 4 
79,366 148.262 


34.929 


25.921 


NORMALTZEO X OISPLACEMENTS 


ANG 


NORMALIZEO VY OISPLACEMENTS 


SINGH 


LUMPEO PARAMETER MEFHOO 


4) 


MOMENTO IN 


yi CORR RUUN NUMBER DF SECTIONS SPECIFIC WT{L8/CU IN) AREA(SO IN) 


x COR CRUE 


30.300909 
625.00000 
625.30090 


4.62000 
29.40000 
29 .40000 


0.28300 
0.28300 
0.28300 
0.28700 
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13.00000 
13.00900 
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NATUFAL FREQUENCIES(RAO/SEC) 


5 
269.188 


25.365 75.162 147.574 


15.285 
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aN 


NORMALIZEO Y OISPLACEMENTS 


SNOW 
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LUMPED PARAMETER METHOO 


MOMENTEIN 4) 


ARPEAISO TN) 


IN) 


SPECTEIC wrtLA/Cis 


NUMBER OF SECTIONS 


Y Be 


x COBROCING 


P?,50090 
625.00000 
4625.90000 


5.88000 
29.40000 
29.40000 


0.28300 
0.28309 
0.28300 
9.28300 


6.00000 
12.00000 
12.00000 
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120.00000 
192.90000 
120.00000 


0.0 
360.00000 


720.00000 


720.00000 
CLAMPEO CLAMPED ENO CONOITION 
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4.62000 
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ROTATORY INERTIA IS NOT INCLUDED 


NATURAL FREQUENCTES(RAD/SEC) 


& 5 
79.284 147.905 210.359 
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25.916 


NORMALIZEO X OISPLACEMENTS 


SNM 


NORMALIZEO Y OETSPLACEMENTS 


SNM TY 


LUMPEG PARAMETER METHOD 


4) 


MIMENTE TA 


Yv a NUMBER OF SECTIONS SPECIFIC wtIts/Cl) IN) ARFACSO IN) 


x COBERCING 
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5.8AN00 
29.40000 
29.40909 


4.00000 0.28300 


13.00000 
13.00000 


120.00000 
192.00000 
120.00000 


0.0 
360.00000 


0.29300 
0.28300 
0.28300 


720.00000 
720.00000 


CLAMPEO CLAMPED ENO 


30.30090 


4.62000 


4.00000 


0.0 


CCNOITION 


RCTATORY INERTIA TS NOT INCLUDED 


NATURAL FREQUENCIES(RAO/SEC) 


Ss 
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FINITE ELEMENT PROGRAM 
General 

Program FINITE is a double precision IBM FORTRAN IV 
language’ digital computer program designed to find the 
principal vibrational modes of straight sided plane frames 
of uniform cross section. The program is written for the 
IBM 360/67 facility installed at the Naval Postgraduate 
School. The program uses a finite-element approach to 
calculate the in-plane undamped natural frequencies and 
normalized displacements of structures. The program neg- 
lects axial and shear deformation, but rotatory inertia 
may be included if desired. 

The program is dimensioned to handle a structure with 
a maximum of 20 members, and a maximum total of 43 elements. 
The program with this dimensioning requires total storage 
of 280K bytes. The properties of each straight member must 
be uniform, but may vary from member to member. 

The more severely constrained end of the structure is 
considered the left end, and must be pinned or fixed. The 
right end may be fixed, pinned, simply supported, or free. 
The program cannot handle any intermediate supports or ex- 
ternally applied forces. 

Program Structure 
The program consists of a main body and four subrou- 


tines. The program requires the use of double precision 


dg 


Sine, cosine, arctangent, absolute value, and square root 
routines. If they are not included in the FORTRAN library 
they must be provided by other means. 

The sections of the main body and their functions are 
explained by comment cards in the program listing. There 
is also a large leading section giving input format infor- 
mation. 

The subroutine ELEM generates the elemental stiffness 
and mass matrices for use in the main program. 

The subroutine CORNER performs the corner transforma- 
tion on the last element to the left of the corner and the 
first element to the right of the corner. 

The subroutines MAPRIN and JACVAT are the same as 
used in the LUMPED program. These subroutines perform the 
same function in the FINITE program as they did in the 
LUMPED program. Once again the rows and columns of the 
mass and stiffness matrices deleted in applying end condi- 
tions are replaced by zeroes. In the mass matrix, each 
diagonal element corresponding to a deleted row and column 
is "repraced wrelmunil ty. 

Input 

The input required by the FINITE program is the same 
as that needed for the LUMPED program with one exception. 
The number of elements into which a member is subdivided 
must be an even number. The input information for the 


LUMPED program is found on page 58 of Appendix B. 
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Output 
The form and type of output for the FINITE program 


is identical to that of the LUMPED program and is discussed 
on page 59 of Appendix B. In order to specify the number 

of frequencies desired, card number 999 in the program list- 
ing is used. The number of frequencies may not exceed twice 


the number of elements used. 
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C-----GENERATION AND COMBINATION OF ELEMENTAL STIFFNESS AND MASS MATRICE 


N 
‘om 
_ © 
NN wv) 
*% € 
% © 
me, 
~ te 
m=) 
oad iy 
fea) = x 
= x 
| ~ -_ 
cam <{ fa 
wd - << 
aja Lu a —m 
n— eee oO ££ 
eed — —_ * 
@ < t= 
(OC) woe ' m=O 
oO ' x~Ow 
0 + T ~O~— 
ce = fa) we NI ened 
tA —_— ot am) al wd PO) 
ox | or ig O < O etl 
QO * a_i CO aw Po WOW 
— - - << ® W~ om Sy, 
eam N Us w= © I< ~ 
ot <{ om mp AON! i= m=O wm e OO 
— | +O. ac <_ ~~ @ -~YT— tom 
eC we on eed Sl on en Od ~~ I mO 
r= UU “> a) + + ~OxKOY 
*¥ +I =—O~—- wZ £ <I OQ >< Lemme wee 
= WO =-O0 Mada bm ond +O - band 7) Km aad CD 
om wee SO “UN @ fb mm LW) <f ® ~~ x rt) 
—<I NE w ek =f FTF pats) ad A) %% ollLlam 
wwe £IOwd SHO AO keh it (ae C”NZ 
ZzaiwxKWNO?W I e + ect UOSE SF a were Zs 
ZW HO om i] = II Wo Sally eM 2£i=— om omYY e@ © 
COC WH PNtt tw oo ee er | ob a 
sa FPF Mt ON Hm tO am = | TH = <Ieded © we ff] wee OY LIL 
LUMA OVI BDZ ems met Ow I LU oe 
ees (FF cat wa? ee we Ze ed ee) et See OW ww ESS lw 
-~ we OTTO Iemma abe Lym WY mm = 
= | ad ed coed oe) oe Sw Cont - es ee 
N Lu OW ZW Bh bt bf Il pred weet fmm wee | oud 


O 
OJ OuUrLoe TOTGCOAQOOO OOK EMD Oestnmu it 
QOEXMNADO we Od Ol DOO Sk Ob ZOQDN SO IO OO tO 


N - © oO oO x 
~~ et wt o=t 


82 


oe | 

o 

= 
OJ a 
0 — Ww 

V— 
= ” Gn 
bn + o« — 

a 

(ae) I= ja © 

+ x< e 
Fi = = cs 
uw = + we «=lCanN 
\ & = | 
CO mm XS OC 
= » ae) rm 
° =e? 8+ & 
oO come) eo) 
e —— om © 
— 
ad lim WM ol 
° i a i € 
= - = el om 
~ + OD em 
+ Fhe ee 
x © Nn = w+} 
~ loan ~ z2Oxy™w 
ad “A NO~x 
< O <I ~~ 
= Oo ~— tirmtg 
= aoe Hee a2 
<I -~ e¢NnZ2as 
So) aa I Dk OO KT 
—w eH eh eo oh D+ NO 
MAO RHONA & mL) 
DMO + od I il i “oem = 0 I} 
<I ON er “we > “| PH YY <T 
eo = wa MwYOwW 


we ee OD OLL TD Il erm LJ 

it bod weer? (9 
wO LOOOrST~VMTUUS 
dad CP) cad tod (2) C2) C1) IN PD PS eet es 


NM oO eS 
00 Se) 


TRANSFORMATIONS 


C~-~--CORNER 


C----- 
C----=- 


K+1) * AR(KK41) * 0.3002) 


| 
— om > eed 
= = hat = =— = 
m ee & = (v) = ee 
Come) oot roms f= + = ond bot 
e eo tes? (/} om rd . awe 
—_ mb > = a pe b= 
— om om (/}/) + & a wv om /)Y) 
be ed bom — +m — te 
am (/) e e+ + Comet Y~Ow mS) e e+ + 
om i om bees > OC we ol = 
OS ee + xX ww CS Ce ed de 
=! = = mh =? a FO am > > fo) 
wom KO HY) — + WO el wee MOC HY) 
= => oe Wow > Ee —} bm 
KY $¢¢7 8 mm YW w& XY +40 
_ mms +> —s bt on 
+2 mame + 7+ ~ + 2 mo + 
pm me tems C7) pe |i | No b— et pO) 
mY Yd mm — — -KY ++Nmm™ 
S) on™ LL LL UL. me etx + Yiw ™ iLL 
ee ee ect eS $Y Om mL eel 


amt el pe te er ee o> a ae Te ee 
bt bh n> Y SeONN Re WDE 


Y= @=MAY ww Ow I Se OR = MWY 
UWI) ere So ee eC Le LY er ee 
—_) > > an) N mnt 7) CO I ee > > i || 
Oz i NAN ~ Ox wi tt NAN 
an _—om “Ii DOD Hem WSN -~m 
cI —- her he % aot = tt ime 
— ii ov”) 7~—— — mem ollie |} $n 
YO on™) bt 379 LentHonwn Yo ome) bm 
CL cam tent fee tat ee NZ $e Ym meine ie | eh 


CL) pm LEN) LN LT mt of XY ee Uk & e+ YY Om o& 
ZY be wee be + we Oh | > OY MM ZZ th + 
CeN SRN MANUN™ © Wi we CUS || eM™ Il 
Cl we ff mt me |) ee LL LL JC tet et eet et et 
Sua) rm mee teed de foc seer wer ed J ere “WL EES be we he er 


MU “MULL “un? x bom WAU WI LL UL 
Sie eae OE el | ede hee © Oe Lee) Sa © 
ceed ed er wae mt et a ee ee eee | 


aErOFzrrOor Fr uHnmu aaEEOF = bh O 
QNMNONNNMOOMN™ = OO O=OONNONNMMG 


19] oO un 
od 


83 


eJJ4+J) + ST(I,J) 
+ XM( I,J) 


N+ + Se 
> © “anew lL 
7) =] Met D 


Ore FTM”AD O 
QW Maz 


4 
5 


Co 


C-----APPLICATION OF ENO CONDITIONS 


C----- 


= 
o-_ 
Ss 
eh SP) 
WN 
aN 
mnt C 
= pole 
w= 
» NO 
~ Oo 
~ o 
vd Q-=— 
ad Nie 
~~ QE 
= mL 
o Oo 
aaa et @ 
rd © 
om © @ 
Bod Ne ot OY 
ot NO Oct 
a & Oo =) 
0 ome ee oD 
od aim Oo OF O0 oo 
ad Na OO e © (ae 
0 ~ vnecectO oaQ es 
wl o<{ O0Qd © OO e 
- Ot ~~ -~ nn oe eed 
Ye) N= CO Mt | om © =a 
=! <{ w= e+ +Oi = Uwe 
ee a> N © eomtTmmTz © oYD> Il Il oo 
oO Q = ath FO FOO oY ae pd =! 
~aA COC HY Met &eO mi Hilt © e+ el I] SD em Il il 
- f= w OQ ANMNAD AN YN weOF DT ey YY YD we 
ODO OO AN Let wt ites NS wwe Ft DD ww OY O 
Oo 28 O00 0x2 © IATUUWLE eZTYOUWLD eDoOM~ 
nL eel | coal oo a | | | | a De ds > ie NP a a a a aN 
peed fox) te? ee? Cee? Qoewd) Qewe) Cees? tee? trae 


P~XOIOIJOBLSLODO“LDPOCDEOOwR eH TE ELOKe TE ZOO 
MMOCOMOmMOKOK MOK KZOOMYNNAIN SOMNUNNN OO 


0 &—-DOO CO - ot a 8) PUY ow Oo QD 
At tet ON NON NN O CN wl 


84 


wW*T=1 6% OA 


-----) 
SYQLISANSODIZA wOdd S3adVHS 300QW i we 0-~~-~2 
(OOT S‘wZ *O305d *T *w *4dIT1S) LVADVE WIVD 
JNNILNOD Ee 
(FSTIASTILS = (CISC IAAIAS 
(C1SfP)ASTLS + (FST) A5T4S) & OGS°O = (FFI) 44TLS 
wW*y=er Ce OG 
I+4¢tfd=» 
W*T=I CE OA ; 
VWANOS 3JNYV XIYLVW SSINSSILS 340 SIN3W3T3S IVNOOVIG 3440 JUNS dai RT 
(V4 WZ & CISFIDIAR «x CIDOIS + (TST) 4S5TLS = (19S E)44ILS ZE 
w*t=f Ze Od 
0°O0 = be OEE IPSS 
w*Tt=I] Ze OA 
W*T=1 2€ OG 
(TIOIA *« (TSF )ODIASZ & CFSTIDASAILS + (T8T)WZ = (78 T)WZ TE 
w*t=f TE OA 
O°O = (15 1)W2 
w*T=I Te GA 
W'T=171 TE OA 
(CEIOTISDLYOSQ / TOTPO = (I1)9IS3 OF 
((I)9I3)S98V0 = (1)9T9 
W'T=I O€ OA , 
XTYLV¥w SSINSATLS JHL NO w3780ud JNIVA “iat 
(OOT *935A3 £9135 ST Sw SWZ) LVADVEF WIV) 
SNNT ANOS 62 
(FST)WZ = CIS )WZ 
(CISfIwZ + (fF S1)IWZ) * OGS°O = (FE E)WZ 
wW*s=f 62 OA 
T #77 = 
W'T=] 62 OA 5 
TWNOSA AYV XTYLVW SSVW NI SIN3W373 IVYNOOVIG 4430 3YNS at Jy Vee 
TOT°O = (IS I)w2 82 
O°0 =H‘ hIw7 Le 
0°0 = (ff I)wWZ 
0O°O0 = (14°F )AATAS 
0°0 = (f*I)45dT4S 


35 


EIVEC(I,tL) + EVEC( I,J) * EIG(J) * ZM(JeL) 


ORT(FREQ(J)) 


© 
@ 
OQ 


TI " 
= £ ZW 
am em aC) 
ont J mt _) emt 
tomy mt tl 


domed eet >) toed 


PUT IN ORDER OF ASCENDING FREQUENCY 


TO 


C---==s5 et ENG 


1 


MN 


we 
Li! 
® 
prea} 
Li 
ee 
N 
ap 
. o) em 
WO tts — 
re awd al 
a8) al ~ = 
aan @) ~~) fal 
Orem Jaw onl 
za te Mw 
= Os Ou e 
i! OOm zo amd 
3 (45 ce ot nom 
on ey >} Wie fe 
© e~ Onsto — 
_ 270M) WwW Aes 
= ot J = «Of — ee | a) 
oe wl aCe | w = Wil tt LL 
mae OD et em e Oe eo 
(CS at en U Yoemelanen Lone atc 
OS VW WL dm IN bm twee Le It ~ e111 OW 
=~ WJ wnt (/) eer (1) DW M— "3 OO <t 
m~ J Wed Wm YS eres 
=) aN QO b> = LU m=O MZ I 
me 1 EO Le Sb SOLU LU ee LJ 
“AD ~~ YO ~>>OOW 
LE YO uw YOu UL Lh UD CS es tat pt bt 
et dC) bee a tet Cl) tet at et fe LO SIL LL LU) 
moO Oo oO w 
OO fS wT a8] 


86 


\ 
| 
j 
j 
J 
O 


C-----X AND Y DISPLACEMENTS 


C----=- 


w™M 
eo 
Sa 
cad Ta 
=" we 
—_ LJ LL 
= = (Lm bet 
—_J @ e —Y VY) 
OOr —e Il ji 
oo on No | ote 
TE | | | el | Tee 
a? ¥ *O Je 
am NS VHEH™ 
ett ACA A AN 
steer try eel 
i= Sar 
O~w~w~ O YY O SU wew 
OK }-O37YIO Ke 
a 
i 


WY) 
™~ 
wT = 
Sg x“ 
o 
© Ww x 
er 6 4 
i 
a8 @ OO 
Go| b= Ls 
> 
—_ Cc) el 
wy Ss) LJ 
| 
om) = + 
- N 
€ doug fh co 
© af 
e= e SY ca 
hex UJ oe = 
—we <r ~~ o 
tT ey ey J ~~ 
id ~~ & oe oo C\J t+ we 
CHc) HO eet em OIC) 
m+ +4 ++x~ ew 
> | awiI~ x > 
CE we tee oe ra at DC er? Dad pot 
oat <qI~ I< | >< Ud 
= = = == + 
m2 md + == il ou 
2<at Hei oi aA <1 Fri 
zone Ow © OOo= -—— a 
Ow pt J ee eT pad www eet pd eed pent 


GIN ee AFC De ODIZM4+O+ $+WWt 
WOIDMStHOTOOsT—O SS IYIDD 
eae Z<Ii<l co 0 ee NOt FH SZQI 
AAN ADOOO 1 OOAOANONA == 


wee & |] ee ww b= = P— & ofp |} 
aa Hn Un Whe weezer 

UL UL. seer LL = © ~Ow~~DOX 

ret at OC al tO OUI SUI OOK OOM 
~ Oo MW Oo wr ON 
sf + vt FT Vv 


C----- 


COORDINATE SYSTEM 


C-----PUTTING DISPLACEMENTS IN REAL 


C----- 


oe) 
~l 


wr 
Ve) 
co 
fon 
CO aoe 
ae) wd om 
o J 
ae —_ & 
GC] Sf) at b= 
p= { < -_ 
oO " x< 
® —”) + 
Te) % 
= ~ WO 
e © —Y”) 
~ e¢ mn > 
ie ©} iy aie, ae 
x _x~y= 
a 6 wv I< el 
= © be ee oo 
Mem LW ~~ Wu 
oJ ¥~ LIYwsz2 
oe od — bh oe 
NN ee ~ On eS 
Vu nO ~KrOZM ££ il Il 
Se S OHOUWIN SOW) 
Teel NYU IOr— 
OO SOO0OrFrO0O00 ot! J 
LUN me ~~ yt 
it "ow Wo ot Dad 
eam tL u. © Z2zaiuww 
OO SMM MN OK De OK 
uw oO O 
uN ON wy 


C----- 


C-----NORMALIZING DISPLACEMENTS 


C----- 


We ele a 
= wee wt J OO 
DOw wei IO Z 
NOK FY OOUOMUW 


N 
uN 


ord 


RETURN TO CALLING 
RETURN TO CALLI- 


= ry 

trond ) Q 

CQ ot N @) 

rad ww 

ad WY 
-~ = rod —") 
> (nl = — 
Ci «4 vd ® ‘ome = 
- ¢ I ee] > cd 
Yr > od Ww — BS © 
—_> het e — ce 
mm WU & 
Li e) i" ry 
— = il ® a 
2 = mal << Ww « 
> ~ << 5 ° 
m fC = x oe ae 
ie (Fa) _ <{ a4 
eo cme fal — fo 
N=-D © bk a8) 
LL mS > Sa <I O OO ma 
> | — - =. C— on & 
OO uW -— oo ww 
a<oee Yom > 6 cm <_{ 
< ot ke o-~ 
Zisz 9 bon be <f = 
ol Ip < > - -_ 
<{~- © CO <f > Tw © << 
w O> WO Od oO oO > 
CO & oH) <f -7 oO - On 
TL, = Ot > o e —~ de 
PH = ond Ws GN — =~) 
WO ow & eZ MM =O 0 e o 
IZ eo = omtnm ~ “YY x oOo pee PP on) 
™) <{ = (1) io mad | oe e Fez 0e LZ I~ ort o 

WtN coe ar e 2 On oe © eedaZ 

LU ot ee UL +t -t Us QS ewted ft if ment } eT eo 
z Z2OeL NN aZ2itezt = ttt=— tisovo 
me ONS te OO eam NWSI DR =MIOKiO 
prae dont toot oom (NJ mer PN OAI~w pe lw on @ Opn d ee ond 


DON eZ ie et DZ eS HOS MON ete 

Ou Zt etka Sede i Ke DOCO™OO~— I aeiw pp Tp 
fo An 8 be Tee eee Pa ee Pee ee ae ot > 48 4 Oe a 
CO Qe SE wee et OF fe tO eee ee pee >>0 wow be OL tae 
DPTZm—LKCOSCWYOWWYOOWLOOO =e OOUWY Om 
Ae OmOU COWL seoemeQOOOWWT OOOO 


= aed 
On MF BON at OS Oe © Oo 
NA NO NS oO OOo OO 
aed emt md owt 


88 


UNEQUAL, 


LW 

ad 
uJ 

x < 
AG 

Sg 
uN 

ee 

‘oa 

& 

ey 

a 
mo 

eo 
Om co oO 
= a oe 

eo _ Con a oon 
we al e & 

<WwW ‘S) ao Mm 

= (oo 

NN — ee 

Ow — e o 

a = m Oo wo 

TW od <q fie Cc, = 

L » oo — we 

od wa »< = _—~ 

= al r= e es ey 
2 eee om oY fot oO Mm WM 

Or = > = = me Vw 
eat be am wt q ww - et YW WM 

os © = Wr <I ~ om 2 
~ = € +4 cc <q ef -~_ 
m— OO -~ ~ + -~ GN eo © = 

eO<ct a oO qa a, ae ee N = - 
T<aq+ > + > ed > = +a * > 
a o— +a oO | Q. <I on ' LY a ee 
co ie) it ome, ~*~ {I = Oo an OQemanwnr w at aN % + 
MOU o _™ Cc 7D ~ re £O 2 QUITE oo ~ ® =x WW 
bond CY bo Te = ©) + QI we & FOr OO = NN Om © 1+ 
See we Ot a NOOO eHOEO | & ~~) oO bh b— 
~urt > Yu< . = com OC (A QA =— de -+ = tH + 
dIay~w~e OC MDM re NIN Ic wsI1I kzenm ee DM AMD fa at g\ am om (UY) 
Wit Ie wD FTF at WFNS WMH OOF — FP CS ert 3 il it 
CEA. ON aeeHHOO TC se OD OF7v7OKIO0On~eg OO OC MNOr!I_OO—- 
fag) eceYty | BOW IP met LS Ot = mow YOm~ ss HO WN Own OCO FF 
w> Wi Dam eRe th Or FO TD —MNHo DWN! Aa ANNE H = MOO 


Kime COOZIO FZO Te + Il holhowOw~— {| Dwr NUN feed & OY 
<TD et et fl eee em he WOM eNO COMA IInd eO OMNCKH Ne ee SD 
Zhe © ob GLO- DI ZF TUL ented ¥ ON Ted Ol WmMIOgma™ Hm OY enwnm © & 
CC eZ ww OZ wr LE OO Need OC BQH Mew eo LOON tlw 
OOCKww OU DeLYOW>S IOONH HNTLOUC—L- PHM WH UOMW UN WU WOH tw 
bh 1) <I <I) oe TW HOLL KT OO ON OO OMG dT Yi et QO HAN OE DAA 
= 


© A Tan Co om wt wy bao oOo o Oo © ~ 
© CO OHA Ona ee oad = ~d a x el N MN 
o—~ m~t detest ote est wo am = ~~ x oa od 4 


89 


*A(TROW, IROW)-S¥*(C¥ATI-S¥AJJ) 
) 
) 


MAO 
% & 
eal 
OO 
OO 
alae, 
= - © a 
~ mt bam on O 
(a8) <I —— N -_ 
fa] % Or of = <I 
cal O om >> * ee WY 
o b= + cA ~~ = On ~— pai 
~ % 3 ™ lam Lt N U, 
N OY m= = i+ mo a O e ~ 
a | + { < % % O b= ad oo, A> + — 
=» © ~wt ez MON + % © — *#t MO ‘ 
~ OW HHWNO o& 1 >= mW” WO =a {I wv KH = = nae 
NN Dee NO Il fl od —-— & = Wy itl NSN Z2oON~ (OO - ee base | & 
oon eae || jj mC Rta CO~m od ec em Ow YW we CS —- = 
~~ NO § wer FO POU oe = WJ om NN t- KS DBD - DBD  oheaal 
a MON OO™ MMH WH BOO A A HIND We et OP NN OP NA 
NO UM ClO ee emnmm Il oO OO OC ae COO mt ee 
SOUS OE mw EN at ee OL tet LD LU ff “LI 
22> AOMere ee OO ef eowD ~— & Dann 1QO~-O0 Dae I~ MDI~ 7 


Cri Sntet FIETCOFSIKHV™!OCL wai VlZlrexeO DO + il fl Zod emdn em eZ 
YL NSN tt OO ONO OF mi NI Se rere DP INO, tt) OO MOM IL 
dant peed fm etree FO eee LOT bm ea A rt 1 etd HO ete OD ete OD 1 > 
mer I Qe Oe ZO oO SKM O INKL Ia> QSOS Hero 
ANLOOC I eww Ot Nee LOO Ha HO tH WE HOON HANTOW tf sO Ile Ow lt Ze 
mm QO DAI Ore DDO Rr WLW NOS QO0OMNOOF Oe aor qtWwodaw 


~~ NA fr - UNO © oO nO + NM 

N NA N OS NA Se) NNO OS Psat aa) 

= as as x i NI sa det oot one mel 
gun 


90 


SUBROUTINE ELEM (XL) 
C-----SUBROUTINE GENERATES ELEMENTAL STIFFNFSS AND MASS MATRICES 


C----- 


C----- 


CS N 
% ca 
+ % 
J ed ed HK = 
x— OKC x 
+t UH % % 


a, ome ome (NI CN) ie a 
AN Net % OANA 
COONOC OAH OOO 
DOHAIDMO aA dIONO 
LNAI Tet ex<LAns 
wa elf e€ «€O rt ee 
eT 0OO # #OO 


om & * © 
% FH = HE 
N A + Ait 
— + N NNQOO at 
om —, t+ + MOA OANUO 
une x Se OO WwW * OO 
© etd x~ XK Sa of e ++ 
on % + O+ 
Tw o + + MQ AH eta edt et 
| aN N= ONONQAGCO#NOAQD 
<I << O =) ct ed OAaQC0OO0OM AMO 
= NO OQ = NODOMST eNO et 
- o AON SON AH SFO eM «€ COWMO ee 
OD om wmd @ @€ Oem @ & Om Om 2 60) 8OO e © ~”-~— 
UVW) OIDNAONONNMOO 2 t+ +O baer Jn 
t+ oe — ff ww ow @& + x + - & 
Ly & HHE NH MH HNE OHH + CNN OH band tent 


awF OOOO M— WW wie MA MMQUW UNM — wo 
Ip x< @e@ @ erdtet beet mde OR NIMOMOOD I ow Mw DO > 
wv) OODMDODQOQ IMO 1OWM I NOOwQwww 1 Ol = Ui x< 


a eur LAW 
Jae eH HH NNN WN NH at om me 
Net Ate 


DD et rd MONIT ON FD SVD A ONO FUL) UD SLD) st Dt ee 


mnt oo” _) ee &® @ F&F e& Fe FF F&F & fF © & & HF © & & O&O O & ~ oO 


me dd et md et QUI ONIN CO CAO SFPUB HIOINICQINQIFO CO COS FUVN CS 


Fm O Ohm mR FS Te Rem eRe Rhee eR ETT SET ESTE TST TOORFTWZ 
= QQVGNAN K KMNMNNN NANNY) KK KOK OE OK OK EK On KOK. CI OKs, LL 


=~ A 


oT 


o 
—_— om 
CJ me 
— 
NSN we 
> 2 —_ 
Lee 
ee 
e 
m LL 
ON = 
_ eo 
p= ae 
= Es 
Looe 
r= 
~- w= 
m WL 
0 & 
- tf 
Tew 
Cam wet ST oe 
uu Lh. er 
= o=— 
- el 
oe | = 
Ke OO onw 
& oe oom Lh 
Uy atwS o& 
ao o~m 
O ~Wu 
oO Zur 
asd & 
o oe am 
Cem ee ome, ome 
ZIG 
~ | fe tone tee’ 
= OW sw 
e ei & 
ZIe N 
ef me alt 
IE 
wew wf © 
eee a 
Zoom 
om OC om Lm 
OH ee I 
a OO eth 
<P 
=<Iwy 
Us <T ee Lis om 
Wk NO-e 
2 fuze 


W 
LUMNS. 


RO 
' Ce 
EMENT FOR F).- 


A 
NK 
AT 


OF MATRIX 
ROWS OF * 
EE OATA ST 


“IY 


VOOVOOO 


- © &® &@ & 
~~. SSS 
By SR, TED EG, 
MrOronsy 

e@9@e¢8@ @ @ 6 
LAW NN OO 
ae) el Hl met et pd 
Ow Ou Ou wv 
BOMOONN 

od at emt mnt || 

- f & & &@ & 
mM OK OK OKO OK LL 
NAN At Net O 

ee fe ee @ e() 
LOLA LA UN UNL CL) 
=) Pe) [=F ed ed ey 
fee? wee Gee ee? ee ee 


~ 
NJ 
ral 
» 
NY = an om on om on LL 
ANNAINOOD 
mond md md od od od CS) 
CS Pad det ng nd bond ng (_) 
MaAOCSOONN me 
rim @ 
C) MK MC KO OK OC 


oGT.6) 


LL LU LL Lu 
Oooehe2 


ANNMMOS TFO-OO0O0ON 
FWA! WAL SS OLLI Oe 
CO eet eee eset ee eet tans evel] () Dom fret bed tad grand 
eELTIELEEE «Cte ot 
ONAN AIOIN RH ORAM += 


J=N1L,N2) 


A(T J)» 


wWUdz <= 


NOVO COCO) CO andere eer eres Eh Tl 2a 
Mm AAR SSS eer LL LL UL UL a | ® ~— ww 
pe LL wad + OO eNYK ISN + O 
OmeNAPFINOOO tt tt tt I ZzOow~wr FTF =F 2 
MULLLILULIL IL Ow rd et ond Lh. LL ead tw 
z UQEawmnenm HZ —-HO CO Z2ezZ 
m= NRHA II bm <T UN bp wc 
<{ <{ el | teeta | i Se Pt | | 
med oe met AIA mad — 2 
<I <f wes TTSTAHNCOOCOCaAVNW!S!S 
20! me SULLA rr ZOU OO OA Zeal 
MAIO SL 


a 


Oo oO 
tr WwW 
ol owl 


Oo 
= 


SUBROUTINE CORNER (ANGLE,NSIDE) 


C----- 


ROTA3 (5,5 


Fa) 
a 
uN 
N 
<j + 
ad 
c= 
cr 
e- 
Ww 
& 
uw 
~ 
a2 
_— @ 
OC 
fad Cc 
bom 
— & if 
Pa) om © 
1unery Loe, ~ 
CO eO 
aw ey a uN 
TC o N Nf /) 
i a = o O 
<I xU 4 QO wt ~N b 
a Cy & tf) — 
oe OD eaedtili C sa ete CC 
D men) AA «= ~ Ff QQ oO 
LEN CD LY tee OO wre e & be = © od ot O@® 
He OA ee @ eC) mm Oeee— e © 
uy oO ODOO0OWeH4 wi | LOOCOrN OO 
aur ZH aS © ele | By 
Th xu) nun woo CONN Un Ne Wd 
iW) eC) = 22 Lu 
(am eTOW LY mom wee cI = od Om ee eel oe 1 OL Oe Fe) 


rm LL) ret beet ee LE ed an) www et TRIAL © ew 
MOM Ost me BF Bw OL test HAZY © OO Bw ODdet OO 
To) eC FH Nl mem meet COCO HO OO ae Sf St mts CUA ff |] 


QNAO — awww eZ 0 8 OO 8 Ow ww www) Sew we 
oad za ANANDTAOOO | OOOO ANNNNS NT agi 
AW2 ANTI w - oa {did md 
az2 A antl end anal cond Wow now | sell aed aoe ened eed and bk 
Zee QQOCOCOW LL O oO OO@OCewoco0 


MOUOAAOOC ALONE NE OWN AEA AOOLe 


N == Qo OO DD es 
— ond 


 ' 


“2 — am ae lige) —— 
— we et ~~ od ond 
ed om ° & MICS ~ & 
I< a ez. qq 57> 
be ~~ —o —~ 
Be oN Toa) OO asm 
er <I << foams 4 <I <{ 
bm be 
% # ele: + # OO 
fae a CY 
= <= Se GFR 
maw Bae | * iar’ Heare | Ht % 
- & - 
— fo) Rey Re ad Cam ey OE 
eee? Gomntt (oe) wa ew wun) Dee==? 
es ao = ~ oe 
WI Tag, >< =+ 
ww wee a? See 
+ + Foon Kool ++ NN 
<I <{ <I <{ 
a OO para Fae | ee] 
ee & oY ~ om@ad 
i) a 
we? mee > + — + + 
NO at (OF) 
- b= aud wand QO FF Sse 
QO@ oo e2« OO ° & 
aw tan od Oo CY ac ad bom 
"ou eo eS Wt Hod oe Fz 
COO NX NM x< 
LD) eam me LN LN ut’ UN LD) ame ome LS) om ome LAL LST uw 
eto eel] i ei i ee eww ee ell il ell il 
et © Od ad od ated © Hed & Hed ead cond 
W] ome mee] if mem ff meme ff | em fl ee i | mm 
Ce ee eee eS 
- & eo oY wt 0) et C(O) ~ & - oY 


ket toes Smt et WOO KT KT OT ee Pe Re et ret Pe es ee 
DD ODOn Ss OnswoOoCOboCOLSs On FW 
OeLrOoOOMNxXOAMXYOOCTAACCOONWNKOMKAW 


on] fr Ww Oo Ls 


94 


FINITE ELEMFNT METHOO 


MOMENT IN 4G) 


TN) ARFA(SO TN) 


NUMBER OF SECTIONS SPECTFIC WTt(LA/CY 


y COGROCTN? 


x COR RGUIN® 


30.300909 
425.00000 
675.9000N 


4.62000 
29.40000 
29.40000 


0.28300 
0.28300 
0.28300 


0.28300 


2.00000 
4.00000 
4.00000 
2.00000 


120.00000 
192.00000 
120.00000 


0.0 
360.00000 
720.00000 
720.00000 
CLAMPEQ CLAMPEO ENO CONOITION 


8?.50909 


S$. 8000 


ROTATORY INERTIA [S NOT INCLUOEO 


NATURAL FREQUENCTES(RAO/SEC) 


S 
280.597 


om 


150.444 


3 
34.679 719.647 


1 
25-752 


NORMALIZEO X OITSPLACEMENTS 


SNM FI 


NORMALIZEO Y OITSPLACEMENTS 


mt OE PUY 


FINTTE ELEMENT METHOO 


MOMENTOCTN 4) 


Y ot VN NUMBER OF SECTIONS SPECIFIC wWrt{Le/Cu TN) AREACSO TN) 


xX COBRRCIN) 


30.30000 
425.90000 
&25.3900N0 


4.62000 
29.40000 
29.40000 


0.28300 
0.28300 
0.28300 
9.23300 


2.00000 
4.00000 
4.00000 
2.00000 


120.00000 
192.00000 
120.00000 


0.0 
360.00000 


720.00000 
720 .00000 
PINNEO PTNNEO ENO CONOTTION 


R2.50000 


§.8R000 


0.0 


ROTATORY INERTTA TS NOT INCLUOEO 


NATURAL FREQUENCIES (RAO/SEC) 


2 4 
75.667 149.820 278.1450 


25.364 


15.256 


NORMALT7EO X OTSPLACEMENTS 


NAO fe ao 

NOaese 

~-OM OO 

SVN == Cc 

eoeeee 

oooeo 
He 


ANAT 


NORMALTZEN Y OTSPLACEMENTS 


KAM 


J 


FINITE ELEMENT METHOD 


MAMPENTIIN 4) 


SPECIFIC WHILR/CU IN) ARFA(SO IN} 


NUMBER OF SECTLONS 


Y COOROTIND 


X COBROCIN) 


82.50090 


625.09009 
625 .90990 


5.-A 8000 
29.40000 
29.40000 


0.28300 
0.28300 
0.29300 
0.28300 


2.00000 
2.00000 
2.00000 


120.00900 
192.00000 
120.00000 


0.0 
360 .00000 


720 .00000 
720 .00000 
CLAMPEO CLAMPED ENO CONDITION 


39.30000 


4.62000 


2.00000 


0.0 


ROTATORY INERTIA IS NOT INCLUDED 


NATURAL FREQUENCTESEURAD/SEC} 


5 
311.542 


34.69% 79.974 152.647 


25.761 


NORMALIZEO X DISPLACEMENTS 


ANNs wh 


NORMALIZED Y DISPLACEMENTS 


=NE SUN 


FINITE ELEMENT METHOD 


&?.50000 
625.90900 
625.30909 


MOMENTECIN 4) 


AREA(SO I} 
5.89000 


SPECLFIC WTI{LB/CU IN) 
0.25300 
0.28300 
0.28300 
0.28300 


NUMRER OF SECTIONS 
4.00000 
4.00000 
4.00000 


Y CORRRCINI 


x COBROTING 


29.40000 
29.40090 
4.42000 


120.00000 


0.0 
360.00000 


720.00000 


720.00000 
CLAMPED CLAMPED END CONDITION 


192.00000 
120.00090 


39.30000 


4.00000 


0.0 


ROTATORY INERTIA IS NOT INCLUDED 


NATURAL FREQUENCIES(RAO/SEC) 


5 
280.594 


6 
79.647 150.444 


2 
34.679 


i 
25.752 


NUKMALIZEO X DISPLACEMENTS 


SNM 


NORMALIZED Y OISPLACEMENTS 


SNM Sin 


96 


APPENDIX D 


TRANSFER MATRIX PROGRAM 

General 

Program DYNRES is a double precision IBM FORTRAN IV 
language / digital computer program designed to obtain 
"exact" solutions for the principal modes of straight 
sided plane frames of uniform cross section. These modes 
are used as a standard for comparing the LUMPED and FINITE 
programs. Axial and shear déformation, as well’ as rotatory 
inertia, are neglected. 

The program uses the method of false position” ime 
assist convergence to the frequencies. Convergence ls gob- 
tained when either the absolute value of the frequency de- 


—756 


terminant is less than 10 , or the frequency is bracketed 


such that the interval is less than veo 

The program is dimensioned for use with structures 
having no more than nine members. The program requires a 
total storage of 50K bytes. Each member must have uniform 
properties, but the properties may vary from member to mem- 
ber. 

The more severely constrained end of the structure is 
designated the left end. This end must be either fixed or 
pinned. The right end is allowed any of the four common 
boundary conditions: fixed, pinned, simply supported, or 
free. The structure may not have any intermediate supports 


or externally applied forces. 


a7 


Program SEructure 


The program consists of the main body and two subrou- 
tines. The program requires the use of double precision 
sine, cosine, arctangent, hyperbolic sine, hyperbolic co- 
Sine, and square root routines. If these routines are not 
available in the FORTRAN library, they must be supplied by 
other means. The main body is prefaced with a comment sec- 
tion containing the format and makeup of the input deck. 
The functions of the various sections of the main program 
are explained in comment cards spaced throughout it. 

The subroutine EVAL generates the transfer matrices 
and a lees the frequency determinant after a trial fre- 
quency has been determined by the main program. 

The subroutine MAPRIN is the output subroutine. It 
performs in the same manner as discussed in the LUMPED 
program on page 58 of Appendix B. 

Input al 

The input required by the DYNRES program is basically 
the same as for the LUMPED program (page 58 of Appendix B). 
The DYNRES program requires the coordinates of the midpoint 
of each member as well as the end and corner coordinates. 
No further subdivision of the frame is required. The rest 
of the data is thé same. Format and order of data cards are 
given at the front of the main program. 

Output 
The coordinates of ends, midpoints, and corners are 


automatically printed by the program. The properties 


oe} 


(modulus of elasticity, specific weight, area, moment of 
inertia) between each set of coordinates are also printed. 
The rest of the output is identical with that of the 
LUMPED program. The number of frequencies solved for is 


specified by a control card. Its makeup is shown in the 


comment section preceding the program. 
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5 0:0 0.7538 0.0479 020435 0.0352 -0.0149 
NORMALIZEO Y OISPLACEMENTS 
l 2 3 4 5 6 
1 0.0 -0:0000 -0.0000 0.3318 0.5511 0.4183 
2 0.0 0.0000 0.0000 0.3911 0.3775 0.1036 
3 0:0 -0:0000 -0-.0000 0.6279 0.0674  -0.6670 
4 0:0 -0.0000 0.0000 0:4770 -0.6017 0.5262 
5 0:0 -0.0000 -0.0000 0.0219 0.0633 <-0.1874 
TRANSFER MATRIX METHOD 
x COGRRCINI Y COOROLIN) MMOULUS OF ELASTICITY(PSI) 
: 30000000.00000 
0.0 60.00000 
30000000.00000 
0.0 120.00000 
30000000.00000 
180.00000 156200000 
" . 30000000.00000 
360.00000 192.00000 
30000000.00000 
540.00000 156200000 
-  30000000.00000 
720.00000 120.00000 
; 30000000. 00000 
720.00000 60.00000 
30000000. 00000 
720200000 iO. 0 
PINNEO PINNEO ENO CONDITION ~ 
ROTATORY INERTIA IS NOT INCLUDEO 
NATURAL FREQUENCIES(RAO/SEC) 
1 2 3 4 
1 15.257 25.365 > 75.652 149.662 
NORMALIZEO X OISPLACEMENTS | 
1 2 3 4 5 6 
I 0.0 0.2982 0.4299 0.4152 0.3940 0.3957 
2 0.0 0.0689 0:0054 °°-0.0927 -0.1257 -0.1777 
3 0:0 0:2040 0:1145 . -0.0153 0.1233 -0.0020 
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3 0:0 . -0:0000 -0:0000 *--0.6493 -0.04638  -0.6702 
4 0:0 -0:0000 0.0000 0.4743 -0.6055 0.5130 
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